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ABSTRACT 


Extensive  measurements  of  thermal  contact  resistance  were  conducted  in  the  following 
study  areas,  where  10  previous  systematic  measurements  have  been  performed:  (1)  Inter¬ 
faces  consisting  of  brittle,  nonmetallic  components,  and  (2)  radiative  transport  across  the 
Interface  at  high  temperatures. 

An  apparatus  for  performing  the  measurements  to  elevated  temperatures  was  designed, 
constructed,  and  performance  tested.  The  equipment,  operated  in  vacuum  environments, 
utilized  tantalum  heaters  to  develop  heat  flow  down  a  specimen  column  which  waS' water 
cooled  at  its  base.  The  contact  resistance  specimen  column  consisted  of  two  high  purity 
beryllium  oxide  cylinders,  1  inch  in  diameter  and  1.5  Inches  in  length.  The  successful 
performance  of  the  apparatus  was  demonstrated  by  measuring  the  thetpial  conductivity  of 
beryllium  oxide  and  tantulum  as  a  function  of  temperature.  The  beryllium  o&de  and  tantalum 
conductivity  results  exhibited  a  precision  of  at  least  ±2%  while  the  accuracy  of  the  data  was 
estimated  to  be  ±3%  by  comparison  with  literature  values. 

Analytical  predictions  of  heat  losses  from  the  specimen  column  were  made  through 
solution  of  the  boundary-value  problem  for  the  temperature  field  in  the  Insulation  segment 
surrounding  this  column.  The  results  compared  favorably  with  experimental  heat  loss 
measurements. 


The  contact  resistance  of  beryllium  oxide  interfaces  with  two  levels  of  surface  roughness 
were  measured:  (1)  a  sandblasted  and  partially  polished  combination  exhibiting  a  nominal 
10  /x  in.  CLA  random-lay  finish,  and  (2)  an  interface  with  sandblasted  surfaces  exhibiting 
a  nominal  40  ^.in.  CLA  random-lay  finish.  Quantitative  surface  roughness  and  waviness 
analyses  were  performed  with  a  Talysurf  Profilometer. 


Contact  resistance  as  a  function  of  mean  interface  temperature  was  measured  for  the 
above  mentioned  combinations  from  about  200°F  to  1600°F  at  5  deadweight  compressive 
loads  ranging  from  5  to  500  psl.  The  reciprocal  of  the  contact  resistance,  the  thermal 
contact  conductance,  exhibited  the  following  trends: 

a.  The  values  for  the  smoother  10  /a  in.  interface  were  higher  by  a  factor  of  3  to  4 
than  those  for  the  40^.  in.  interface  at  comparable  compressive  load  levels. 
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b.  The  values  Increased  rapidly  with  compressive  load. 

c.  Radiative  heat  transport  across  the  interfaces  was  important  at  high  temperatures, 
particularly  at  low  compressive  load  levels.  This  contribution  increases  by  a  factor  of  up 
to  20  in  the  temperature  range  from  200  to  1600aF  and  was  a  major  contributor  to  the  net 
interfacial  heat  transport. 

d.  The  compressive  load  on  the  interface  serves  to  reduce  the  relative  contributions 
of  radiative  transport  indirectly  by  preventing  the  generation  of  large  interfacial  temperature 
drops  at  higher  temperatures. 

A  combined  roughness-waviness  interfacial  heat  transfer  model  was  formulated  which 
successfully  correlated  the  experimental  measurements  within  10%,  including  the  contri¬ 
butions  of  high  temperature  radiative  transport  across  the  contact.  The  thermal  resistance 
induced  by  surface  waviness  were  in  series  with  those  induced  by  the  roughness,  the  com¬ 
bination  being  in  parallel  with  the  thermal  resistance  across  the  radiation  gap  regions. 

The  nonlinear  character  of  the  radiative  transport  effects  were  successfully  incorporated 
in  a  finite-element  temperature  field  analysis  that  exhibited  stable  convergence  over  wide 
ranges  of  radiation  heat  transfer  coefficient. 


The  relative  contributions  of  radiative  and  conductive  transport  across  the  interface 
and  their  combined  influence  on  the  temperature  field  distortion  were  characterized  by  a 
dimensionless  interfacial  heat  transport  modulus,  =  h  k.  This  group  was  proportional 
to  the  4.4  power  of  the  mean  interface  temperature,  T  ,  expressed  in  °R;  a  4.0  dependence 
was  predicted  theoretically  for  the  beryllium  oxide  interface. 


A  simplified  analytical  procedure  for  predicting  the  temperature  dependency  of  the 
thermal  contact  conductance  was  also  developed  and  applied  successfully  giving  essentially 
the  same  results  as  the  more  detailed  combined  model  analysis.  This  simplified  procedure 
was  partially  empii^al  in  nature,  however. 


iv 


AFML-TR-69-1 


TABLE  OF  CONTENTS 

SECTION  '  PAGE 

I  INTRODUCTION  1 

II  APPARATUS  CONSTRUCTION  AND  QUALIFICATION  TESTING  5 

A.  Equipment  Development  5 

B.  Design  Considerations  13 

1.  Design  Philosophy  13 

2.  Instrumentation  13 

3.  Heat  Loss  Analysis  17 

C.  Apparatus  Performance  24 

1.  Thermal  Conductivity  Specimens  24 

2.  Thermal  Conductivity  Measurements  on  BeO  27 

3.  Thermal  Conductivity  Measurements  on  Tantalum  30 

4.  Heat  Loss  Measurements  33 

III  CONTACT  RESISTANCE  RESULTS  37 

A.  Contact  Surface  Characterization  37 

B.  Measurement  Procedure  42 

C.  Graphical  Presentation  and  Qualitative  Interpretation  of  Results  45 

IV  DISCUSSION  OF  RESULTS  60 

A.  Formulation  of  the  Contact  Model  60 

1.  General  Contact  Characteristics  60 

2.  Surface  Roughness  62 

3.  Surface  Waviness  64 

4.  The  Combined  Model  73 

B.  Comparison  of  Analytically  Predicted  Results  with  Experimental 

Measurements  76 

1.  Input  Parameters  76 


v 


AFML-TR-69-1 


TABLE  OF  CONTENTS  (Contd) 

SECTION  PAGE 

2.  Analytically  Developed  Contact  Conductances  78 

C.  A  Simplified  Analytical  Prediction  Procedure  87 

V  CONCLUSIONS  AND  RECOMMENDATIONS  90 

A.  Conclusions  90 

B.  Recommendations  92 

APPENDIX  I.  HEAT  LOSS  ANALYSIS  97 

A.  Solution  Considering  Boundary  Condition  4  97 

B.  Solution  Considering  Boundary  Conditions  1  and  2  102 

C.  Evaluation  of  Fourier  Coefficient  Integrals  108 

APPENDIX  II.  FINITE  ELEMENT  ANALYSIS  118 

A.  Nodal  Equation  Derivations  118 

B.  Numerical  Evaluation  125 

APPENDIX  III.  PHYSICAL  PROPERTY  DATA  127 

A.  Thermal  Conductivity  of  the  Alumina  Bubble  Insulation  1 27 

B.  Beryllium  Oxide  Emittance  Data  1 29 

APPENDIX  IV.  TABULATED  EXPERIMENTAL  DATA  134 

A.  Beryllium  Oxide  Thermal  Conductivity  Data  134 

B.  Tantalum  Thermal  Conductivity  Data  135 

C.  Thermal  Contact  Resistance  Data  135 

REFERENCES  146-148 


vi 


AFML-TR-69-1 


LIST  OF  ILLUSTRATIONS 

FIGURE  PAGE 

1.  Schematic  of  Cut-Bar  Comparative  Conductance  Apparatus  6 

2.  Contact  Resistance  Apparatus  Cross-Section  8 

3.  Assembled  Components  10 

4.  Overall  View  of  the  Control  -  Instrumentation  Assembly  11 

5.  Assembled  Test  Cell  Including  Load  Column  12 

6.  Perturbation  of  Heat  Flow  From  a  Nonconducting  Cylinder  in  an 

Isotropic  Media  15 

7.  Temperature  Field  Perturbation  as  a  Function  of  Specimen 

Conductivity  17 

8.  Insulation  Segment  and  Boundary  Conditions  18 

9.  Typical  Temperature  Distributions  During  BeO  Thermal 

Conductivity  Measurements  28 

10.  Thermal  Conductivity  of  Beryllium  Oxide  as  a  Function  of 

Temperature  29 

11.  Thermal  Conductivity  of  Beryllium  Oxide  Plotted  as  a  Function 

of  (1/T)  31 

12.  Thermal  Conductivity  of  Tantalum  as  a  Function  of  Temperature  32 

13.  Surface  Characterization  Equipment  38 

14.  Surface  Roughness  Profile  Traces  39 

15.  Surface  Waviness  Profile  Trace  40 

16.  Solid  Area  as  a  Function  cf  Vertical  Displacement  41 

17.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

40^i.in./40 /j. in.  Interface,  5. psia  Compressive  Load  46 

18.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

40  fj. in./40  /x in.  Interface,  65.2  psia  Compressive  Load  47 

19.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

40  ^xin./40  /x in.  Interface,  149.3  psia  Compressive  Load  48 

20.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

40  fitn./40  fiin.  Interface,  368.6  psia  Compressive  Load  49 

vii 


AFML-TR-69-1 


LIST  OF  ILLUSTRATIONS  (Contd) 


FIGURE 

21.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

40  fxicu/iO  pin.  Interface,  511.5  psia  Compressive  Load 

22.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

10  pinJlQ  ^.In.  Interface,  4.43  psia  Compressive  Load 

23.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

10  ^i.in./10  ^tin.  Interface,  64.7  psia  Compressive  Load 

24.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

10  /itn./10  ^tln.  Interface,  148.7  psia  Compressive  Load 

25.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

10  p In./l0^itn.  Interface,  368.5  psia  Compressive  Load 

26.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

10  plnJlQ  pin.  Interface,  512.0  psia  Compressive  Load 

27.  Vacuum  and  Inert  Gas  Environment  Contact  Resistance  Results: 

40  p in./40  pin.  Interface 

28.  Vacuum  and  Inert  Gas  Environment  Contact  Resistance  Results: 

10  f iln./10  pin.  Interface 

29.  Flux  Field  Distribution  at  an  Interface 

30.  Idealized  Contact  Element  and  Associated  Heat  Flux  -  Temperature 
Field 

31.  Contact  Element  Quadrant 

32.  Finite  Element  Nodal  Point  Definitions 

33.  Isotherm  Distribution  in  the  Contact  Element  Quadrant  with  Normal 
Radiation  Heat  Transfer  Coefficient 

34.  Isotherm  Distribution  in  the  Contact  Element  Quadrant  with  Increased 
Radiation  Heat  Transfer  Coefficient 

35.  Combined  Model  Resistance  Analog 

3S.  Experimental- Analytical  Contact  Conductance  Comparisons:  10  ^.tn. 

Interface 

37.  Experimental-Analytical  Contact  Conductance  Comparisons:  40  ^.in. 

Interface 

38.  Interfaciai  Heat  Transport  Modulus  for  the  10  and  40  pin.  Interfaces 


T9 


AFML-TR-69-1 


LIST  OF  ILLUSTRATIONS  (Contd) 


FIGURE 

39.  Experimental-Analytical  Contact  Conductance  Comparisons: 
Simplified  Analysis  -  10  fiin.  Interface 

40.  Experimental-Analytical  Contact  Conductance  Comparisons: 
Simplified  Analysis  -  40  /a in.  Interface 

41.  Thermal  Conductivity  of  Aluminum  Oxide  Powder  Insulations 

42.  Spectral  Normal  Emittance  of  Beryllium  Oxide  in  the  Infrared 

43.  Total  Normal  Emittance  of  Beryllium  Oxide  as  a  Function  of 
Temperature 


PAGE 

88 

89 

128 

130 

132 


ix 


AFML-TR-69-1 


LIST  OF  TABLES 

TABLE  PAGE 

I  Positive  Roots  of  u  (A'njr  ;  r. )  =  0  When  r./r  =  6.0  24 

o  c  o  l  10 

II  Composition  of  Beryllium  Oxide  Test  Specimens  25 

III  Composition  of  Tantalum  Test  Specimens  26 

IV  Experimental  -  Analytical  Heat  Exchange  Comparison  34 

V  Contact  Resistance  Error  Analysis  Results  44 

VI  Nodal  Equation  Summary  68 

VII  Effective  Gap  Width  Comparisons  74 

VIH  Surface  Geometry  Parameters  77 

tic  Roughness-Waviness  Contact  Resistance  Results  81 

X  Interfacial  Heat  Transport  Modulus  Results,  Hj  82 

XI  BeO  Thermal  Conductivity  Guarding  Data  1 1 5 

XII  Adjusted  Temperature  Distribution  Data  116 

XIII  “First  Guess”  Temperature  Distribution  125 

XIV  Nodal  Equations  for  Computer  Evaluation  126 

XV  Beryllium  Oxide  Emittance  Specimen  Characteristics  133 

XVI  BeO  Thermal  Conductivity  Results  134 

XVII  Tantalum  Thermal  Conductivity  Results  135 

XVm  Thermal  Contact  Conductance  Results  -  40 ^xin.  /40/iin. 

Interface:  5.  02  psia  CompresBive  Load  136 

XIX  Thermal  Contact  Conductance  Results  -  40  ^.in.  /40^,in. 

Interface:  65.  2  psia  Compressive  Load  137 

XX  Thermal  Contact  Conductance  Results  -  40yu,in. /40^.in. 

Interface:  149. 3  psi  Compressive  Load  138 

XXI  Thermal  Contact  Conductance  Results  -  40  ^tin. /40  ^tin. 

Interface:  368. 6  psi  Compressive  Load  139 

XXII  Thermal  Contact  Conductance  Results  -  40  ^.in,  /40^in. 

Interface-  511.  5  psi  Compressive  Load  140 


x 


AFML-TR-69-1 


LIST  OF  TABLES  (CONTD) 


TABLE 

PAGE 

XXIII 

Thermal  Contact  Conductance  Results  -  10  pin.  /10  ^tin. 

Interface:  4.43  psla  Compressive  Load 

141 

XXIV 

Thermal  Contact  Conductance  Results  -  10  ^tln. /lO^tin. 

Interface:  64.7  psl  Compressive  Load 

142 

XXV 

Thermal  Contact  Conductance  Results  -  lO^tin. /10 /i.in. 

Interface:  148.7  psla  Compressive  Load 

143 

XXVI 

Thermal  Contact  Conductance  Results  -  10  ^.ln. /lO^tln. 

Interface:  368.  5  psia  Compressive  Load 

144 

XXVII 

Thermal  Contact  Conductance  Results  -  lO^tin. /lO^tin. 

Interface:  512.  0  psia  Compressive  Lead 

145 

AFML-TR-69-1 


NOMENCLATURE 


a 

A 

b 

C 
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G(rt;  r^) 


HI 
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I  ,K 
o*  o 


J  ,  Y 
o*  o 

k 

1 

L 

D 

q 

Q 

r 

R 
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<K^n:  r,i  rj) 


solid  contact  area  radius,  in. 

2 

cross-sectional  area,  ft 
adiabatic  contact  element  radius,  in. 
heat  capacity  at  constant  volume,  BTU/lb-'F 
average  particle  diameter,  microns 
geometrical  factor 

2 

blackbody  emissive  power,  BTU/hr-ft 
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constriction  alleviation  factor  defined  in  Eq.  2P 

cylindrical  function  defined  in  Eq.  9 

o 
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interfacial  heat  transport  modulus 
finite  element  generalized  indicies 

zero  order  modified  Bessel  functions  of  the  first  and  second  kind, 
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zero  order  Bessel  functions  of  the  first  and  second  kind,  respectively 
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resistance  product  defined  in  Eq.  32,  (hr-*F/hr) 
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eigenvalue  parameter 

eigenvalue  parameter 
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SECTION  I 

INTRODUCTION 


When  two  surfaces  are  brought  together  to  form  an  Interface  the  true  soltd-to-soltd 
contact  area  between  them  is  generally  a  small  fraction  of  the  apparent  area  over  which 
they  meet.  This  direct  contact  area  may  be  less  than  1%  of  the  total  and  rarely  exceeds 
about  10%  unless  bonding  agents  are  introduced.  It  is  intuitive  that  when  an  Interface  is 
encountered  in  a  heat  flow  field  there  will  be  a  resistance  to  heat  transfer  across  it.  As  a 
result  a  temperature  discontinuity  will  appear  across  the  contact;  this  temperature  drop 
divided  by  the  heat  flux  flowing  across  the  interface  is  termed  the  thermal  contact  resistance. 

Early  work  in  the  contact  resistance  field  was  centered  on  the  study  of  electrical  energy 
flow  in  various  types  of  electrical  motor  and  switch  applications.  More  recently  a  great 
amount  of  activity  has  been  underway  in  studying  the  flow  of  heat  from  metal-clad  nuclear 
fuel  elements.  Contact  resistance  is  very  important  here  because  of  the  large  heat  fluxes 
crossing  the  fuel  element-cladding  interface.  Other  current  applications  where  Interfacial 
heat  transfer  is  important  include  high  temperature  aircraft  structural  design,  dissipation 
of  heat  from  space  vehicle  systems,  heat  exchanger  design,  and  the  design  of  turbine  powered 
engines  and  solid  and  liquid  fueled  rockets.  Recent  Soviet  activity  in  the  area  has  emphasized 
these  latter  applications. 

An  example  of  the  practical  importance  of  joint  thermal  conductance  in  aerospace 
structural  component  performance  is  illustrated  in  the  work  of  Griffith,  et  al.  (Reference  1) 
who  studied  the  structural  response  of  skin-stiffener  panels  heated  on  their  external  surfaces. 
These  constructions  are  typical  of  aircraft  and  space  vehicle  structural  units.  The  presence 
of  joints  greatly  altered  temperature  distributions  and  increased  thermal  stresses  in  the 
panels.  In  a  Mach  2  jet  tunnel  test  of  a  wing  section  it  was  observed  that  flutter  and  dynamic 
failure  resulted  from  poor  thermal  contact  in  the  Joint. 

In  higher  temperature  aerospace  applications  such  as  rocket  nozzles  and  high  perfor¬ 
mance  reentry  vehicles,  many  of  the  materials  used  are  oxidation  resistant  nonmetallic 
refractories.  In  these  systems  interfaces  are  often  encountered  which  may  constst  of  metallic/ 
nonmetallic  or  nonmetallic/nonmetalllc  components.  During  system  operation  it  is  generally 
found  that  heat  is  transferred  across  these  interfaces  at  high  temperature  levels.  A  knowledge 
of  thermal  contact  resistance  ranges  for  such  high  temperature  interfaces  is  clearly  neces¬ 
sary.  This  is  particularly  true  for  the  case  of  nonmetallic  interfaces  where  the  materials 
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are  of  a  brittle  nature.  Accurate  thermal  stress  analysis  if  vital  in  design  with  brittle 
materials  and  it  is  only  with  accurate  estimates  of  the  temperature  fields  in  complex  struc¬ 
tures  that  these  thermal  stress  calculations  can  be  made.  Examples  of  specific  engineering 
applications  in  this  area  include  the  study  of  interface  contact  coefficients  in  the  Minuteman 
solid  rocket  nozzle  (Reference  2)  and  the  contemplated  use  of  nonmetallic  thermal  contact 
resistance  barriers  as  a  form  of  insulative  liner  for  high  pressure  thrust  chambers. 


In  a  typical  solid  rocket  the  upper  stages  utilize  a  flame  barrier  backed  up  by  a  heat 
sink  of  graphite.  In  connection  with  this  application  Gilchrist  (Reference  la)  studied  heat 
transfer  across  contacts  consisting  of  graphite,  tungsten,  and  molybdenum,  including  various 
combinations  of  these  materials.  Interface  temperatures  from  1000  to  2000nF  were  studied 
in  this  investigation,  together  with  high  temperature  bonding  and  brazing  schemes.  In  general, 
the  results  Indicated  that  interfacial  temperature  drops  from  about  100  to  1000°F  were  likely 
at  the  heat  flux  levels  encountered  in  such  high  performance  rocket  systems.  The  possible 
thermal  stress  and  reliability  uncertaities  can  become  considerable  in  such  instances. 


The  author  has  prepared  Reference  1  with  the  objective  of  reviewing  and  summarizing 
the  literature  in  the  thermal  contact  resistance  field  with  particular  attention  to  high  tem¬ 
perature  results  and  data  on  nonmetallic  interfaces.  The  volume  also  reviews  in  some 
detail  the  available  analytical  approaches  used  for  both  single  and  multiple  contact  interfaces 
and  the  more  successful  experimental-analytical  correlations.  The  correlations  developed 
as  well  as  the  literature  results  which  were  of  direct  value  will  be  introduced  and  discussed 
in  later  portions  of  this  report. 


Several  general  observations  relevant  to  high  temperature  applications  can  be  made  as 
a  result  of  the  literature  review  (which  has  been  continually  updated  to  the  time  of  this 
writing): 


a.  The  overwhelming  majority  of  thermal  contact  resistance  studies  have  been 
concerned  with  metal 1  i  c/me  tall  ic  interfaces. 

b.  The  independent  variable  most,  often  employed  is  the  compressive  load  on  the 
Interface. 

c.  In  those  few  cases  where  contact  resistance  is  measured  as  a  function  of  tem¬ 
perature,  the  maximum  temperature  ranges  between  500-700'F. 
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Based  on  the  practical  importance  of  contact  resistance  data  for  high  temperature 
aerospace  applications  and  the  paucity  of  such  information  as  reflected  in  the  above  obser¬ 
vations,  a  program  was  initiated  to  study  the  contact  resistance  characteristics  of  non- 
metallic  interfaces  extending  to  elevated  temperatures.  To  the  best  of  the  author's  knowledge 
no  previous  systematic  study  of  heat  transfer  across  nonmetallic  interfaces  has  been  made 
nor  have  there  been  any  previous  investigations  conducted  in  the  range  though  necessary  to 
account  for  significant  radiative  contributions  to  the  overall  interfacial  heat  transfer. 


The  selection  of  the  nonmetallic  materials  for  study  was  based  primarily  on  their 
aerospace  importance  and  on  experimental  requirements  for  high  accuracy  contact  resistance 
measurements.  The  oxide  ceramics  were  selected  initially  for  consideration  because  of 
their  broad  technological  usage.  From  an  experimental  point  of  view,  it  was  found  that  the 
level  of  thermal  conductivity  was  one  of  the  most  important  parameters,  high  conductivity 
being  desirable  for  highest  measurement  accuracy.  This  criteria  lead  immediately  to  the 
selection  of  beryllium  oxide  which  has  the  highest  conductivity  of  the  oxide  ceramics. 

The  material  has  the  additional  attribute  of  being  reproducibly  available  as  high  purity 
stock  from  commercial  vendors;  its  intrinsic  toxicity  is  minimal  in  sintered  form  and,  with 
reasonable  precautions  in  handling,  presents  no  problems.  Further,  the  temperature  de¬ 
pendence  of  the  thermal  conductivity  obeys  the  classical  hyperbolic  functionality  expected 
theoretically  for  pure  phonon  conducting  dielectrics.  Because  the  berry  Ilium  oxide  specimens 
were  utilized  as  heat  meters  in  the  thermal  contact  resistance  measurements,  a  precise 
knowledge  of  the  conductivity  was  necessary. 


The  major  variables  in  the  thermal  contact  resistance  measurement  program  on  the 
beryllium  oxide  interfaces  were:  (1)  the  surface  finish,  (2)  the  mean  interfacial  tempera¬ 
ture,  and  (3)  the  compressive  load. 


Section  II  describes  the  equipment  constructed  for  this  investigation  and  summarizes 
the  results  of  the  experimental  checkout  of  the  equipment,  which  included  thermal  conductivity 
measurements  on  beryllium  oxide  and  tantalum.  Comparison  of  these  results  with  literature 
values  were  made  in  assessing  equipment  performance.  Analytical  equations  were  also 
developed  in  this  section  for  predicting  heat  losses  in  the  equipment;  full  derivations  of 
these  equations  are  given  in  Appendix  I.  The  predicted  heat  losses  were  then  compared  with 
experimentally  measured  values  a/;  a  further  evaluation  of  apparatus  performance. 
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The  next  major  section  of  the  report  (Section  ID)  presents  the  experimental  thermal 
contact  resistance  results  along  with  important  surface  characterization  information. 

Following  this  is  the  discussion  of  the  results,  Section  IV,  which  includes  the  develop¬ 
ment  of  two  analytical  procedures  for  the  prediction  of  thermal  contact  resistance.  Again, 
the  details  of  the  analytical  derivations  were  relegated  to  an  appendix.  Detailed  comparison 
was  made  between  these  predicted  contact  resistance  curves  and  the  experimentally  obtained 
results. 

Section  V  covers  conclusions  and  recommendations. 
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SECTION  11 

APPARATUS  CONSTRUCTION  AND  QUALIFICATION  TESTING 

A.  EQUIPMENT  DEVELOPMENT 

One  of  the  most  Important  problems  encountered  in  conducting  experimental  contact 
resistance  measurements  is  controlling  the  wariness  of  the  test  interface.  As  amply  demon¬ 
strated  in  the  literature  (Reference  1).  inconsistencies  and  nonreproducibility  of  results  can 
very  often  be  traced  to  this  factor,  particularly  at  low  compressive  load  levels.  In  addition 
to  careful  surface  preparation,  one  of  the  simplest  methods  of  minimizing  this  problem  is 
to  work  with  interfaces  of  minimum  cross-sectional  area.  Another  Important  consideration 
for  thermal  contact  resistance  measurements  is  the  development  of  a  uniform  heat  flux 
field  on  both  sides  of  the  interface  to  facilitate  beat  flux  and  terr  perature  measurements. 
Meeting  these  two  requirements  leads  to  specimen  configurations  with  length  to  diameter 
ratios  greater  than  unity.  Thus,  in  well  designed  contact  resistance  experiments,  cylindrical 
specimens  with  length  to  diameter  ratios  from  1.5  to  3.0  are  widely  used. 

For  such  geometries,  axial  heat  flow  is  highly  advantageous  as  long  as  the  level  of 
specimen  thermal  conductivity  is  fairly  high,  “High”  in  this  sense  lmplioe  relative  to  the 
conductivity  of  the  insulation  used  along  the  lateral  surfaces  of  the  specimens  to  prevent 
radial  heat  losses.  The  equipment  employed  in  this  investigation  is  a  modification  of  the 
“cut  bar”  design  used  widely  in  thermal  conductivity  measurements  on  high  conductivity 
materials  with  cylindrical  specimen  configurations. 

The  basic  design,  shown  schematically  in  Figure  1,  is  quite  simple.  A  composite  bar 
of  circular  cross  section  is  built  up  of  several  specimens  by  stacking  one  upon  another. 
One  end  of  the  column  is  heated,  the  other  cooled  to  produce  a  steady  axial  heat  flux  down 
through  the  composite  bar.  The  column  is  insulated  laterally,  and  guard  heaters  are  employed 
to  reduce  radial  heat  flow.  The  series  of  guard  heaters  are  used  to  match  the  axial  tem¬ 
perature  distribution  at  the  external  surface  of  the  insulation  to  that  down  through  the 
specimens.  If  the  radial  heat  losses  are  small  then  at  steady  state  the  ratio  of  (dT/dz)  in 
each  specimen  is  directly  equal  to  the  ratio  of  the  respective  thermal  conductivities.  In 
comparative  conductivity  determinations ,  heat  meters  with  known  conductivities  are  placed 
above  and  below  the  specimen  on  which  thermal  conductivity  is  to  be  determined  (Figure  1). 
A  nominal  compressive  load  (10  to  100  psi)  is  generally  applied  to  reduce  the  contact  re¬ 
sistance  between  specimens. 


X -THERMOCOUPLES 


Figure  1.  Schematic  of  Cut-Bar  Comparative  Conductance  Apparatus 
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For  contact  resistance  measurements,  two  specimens  of  known  conductivity  are  brought 
together  to  form  the  interface.  The  temperature  gradient  in  the  specimens  away  from  the 
interface  and  the  known  thermal  conductivity  are  sufficient  to  determine  the  heat  flux  crossing 
the  interface.  It  is  obvious  that  the  relative  magnitude  of  lateral  heat  losses  will  be  small 
if  the  specimen  conductivity  is  high.  The  interfacial  temperature  drop  generated  at  the  mutual 
contact  between  the  specimens  is  determined  by  extrapolating  the  steady  state  temperature 
gradients  in  the  specimens  to  the  contact  plane.  The  thermal  contact  resistance  is  then 
defined  as  the  ratio  of  this  interfacial  temperature  drop  to  the  steady  state  axial  heat  flux 
in  the  column.  Provision  is  made  for  application  of  a  variable  dead  weight  compressive  load 
along  the  axis  of  the  specimens. 

The  particular  design  used  in  this  investigation  has  several  distinctive  features  which 
will  be  brought  out  in  discussing  the  apparatus  cross  section  shown  in  Figure  2.  Overall, 
the  apparatus  has  been  designed  for  high  temperature  operation.  The  specimens  (1-inch 
diam.  x  2-inch  long  cylinders)  for  either  thermal  conductivity  or  thermal  contact  resistance 
measurements  are  centered  on  a  water-cooled  copper  pedestal.  The  pedestal  is  machined 
internally  with  channels  both  in  the  upper  face  and  in  the  neck  such  that  water  circulates 
upward  from  the  bottom  of  the  pedestal,  reverses  it  path,  and  returns  downward.  With  this 
arrangement  no  water  connections  are  required  within  the  vacuum  chamber  which  surrounds 
the  entire  unit  and  which  rests  on  the  base  plate. 

The  specimens  are  insulated  laterally  with  low  conductivity  alumina  bubble  insulation 
which  is  contained  within  the  tantalum  cylinder.  The  tantalum  cylinder  has  a  flange  near  its 
upper  extremity  which  supports  the  guard  heater.  The  lower  end  of  this  cylinder  is  supported 
and  cooled  by  the  copper  pedestal.  The  stainless  steel  and  tantalum  radiation  shields  for  the 
guard  heater  rest  on  a  ring  of  transite  insulation.  This  ring  is  supported  from  the  base  plate 
by  an  Inconel  702  support  stand.  The  upper  portion  of  the  Inconel  stand  is  polished,  thus 
serving  as  an  additional  radiation  shield. 

The  main  heater  core  is  supported  by  an  alumina  disk  which  in  turn  rests  on  the  alumina 
bubble  insulation.  This  disk  also  serves  to  support  the  five  tantalum  radiation  shields  for 
the  main  heater.  Because  of  the  high  operating  temperature  of  the  heaters  it  was  calculated 
that  radiation  shields  would  be  more  efficient  thermally  than  other  forms  of  insulation. 

A  dead  weight  compressive  load  is  transferred  to  the  specimen  column  through  a  tantalum 
plug  which  absorbs  heat  from  the  main  heaters  but  floats  independent  of  it.  The  stainless 
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steel  load  plunger  is  thermally  insulated  from  the  tantalum  plug  by  a  section  of  high  com¬ 
pressive  strength  low  thermal  conductivity  zlrconia  (calcta  stabilized). 

Instrumentation  with  platinum/platinum- 10%  rhodium  thermocouples  is  shown  in  Figure 
2.  In  this  particular  setup  the  three  specimens  have  thermocouples  attached  to  their  surfaces. 
The  leads  of  these  surface  couples  are  electrically  Insulated  from  the  specimen  and  wrapped 
circumferentially  about  one-half  turn  along  an  isothermal  plane  to  minimize  heat  leaks  down 
the  leads.  Additional  thermocouples  are  placed  in  radial  holes  drilled  to  the  centerline  of 
each  specimen.  These  couples  are  placed  in  two-bore  ceramic  tubing  which  extends  radially 
from  the  centerline  out  beyond  the  surface  of  the  specimen  into  an  Isothermal  plane;  again, 
this  is  to  prevent  heat  leak  down  the  leads  of  the  couples.  High  temperature  "Astroceram" 
cement  (from  American  Thermocatalytic  Corp. ,  Mineola,  New  York)  was  used  to  anchor  the 
various  couples. 

During  operation,  the  main  heater  power  is  set  to  provide  the  desired  temperature 
distribution  in  the  specimen  column.  The  guard  heater  is  adjusted  so  that  the  temperature 
difference  between  thermocouples  “C”  and  “D”  (Figure  2)  is  zero.  In  some  instances  these 
couples  were  connected  differentially  to  allow  automatic  control  of  the  guard  heater. 

Figure  3  gives  a  closeup  photograph  of  the  assembled  apparatus.  In  Figure  4  is  presented 
an  overall  view  of  the  equipment.  The  water-cooled  vacuum  chamber  is  lowered  onto  the 
base  plate  during  vacuum  or  inert  environment  testing.  The  four  upright  posts  seen  in 
Figure  4  support  the  compressive  loading  assembly.  This  loading  assembly  and  the  vacuum 
chamber  are  shown  in  place  in  Figure  5. 

The  fulcrum  of  the  compression  unit  is  mounted  in  a  heavy  gage  ball  bearing  block. 
The  counterweights,  the  compressive  load  ant  the  main  lever  arm  are  all  loaded  on  knife 
edges.  The  lever  arm  has  a  16  to  1  mechanical  advantage.  The  exact  dimensions  of  the 
lever  arm  makes  it  simple  to  calculate  the  compressive  load  established  by  the  plunger  as 
a  function  of  the  weight  on  the  load  pan.  Experimentally  this  was  checked  with  the  aid  of  a 
proving  ring  to  a  load  level  of  5200  psi.  Counterweights  were  first  added  to  balance  the  arm 
at  zero  compressive  load;  following  this,  the  actual  load  was  measured  as  a  function  of  pan 
weight. 

A  deviation  plot  of  the  results  was  prepared  by  taking  the  percentage  difference  between 
calculated  and  measured  values  at  various  compressive  load  levels  and  then  plotting  these 
values  as  a  function  of  load  level.  Deviation  from  theory  was  within  2%,  while  the  deviation 
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Figure  3.  Assembled  Components 
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Figure  5.  Assembled  Test  Cell  Including  Load  Column 
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from  an  average  curve  through  the  data  was  within  about  1%.  At  load  levels  of  1450  pel  and 
4200  psi,  runs  were  made  in  which  increases  and  then  decreases  in  load  about  these  values 
were  used  to  assess  the  effect  of  friction  within  the  assembly.  Data  points  for  increasing  and 
decreasing  loads  at  these  two  levels  were  within  1%,  and  thus  frictional  effects  appear  to  be 
quite  small. 

B.  DESIGN  CONSIDERATIONS 

1.  Design  Philosopy 

The  Equipment  layout  ]ust  described  assists  in  the  following  discussion  of  the 
rationale  which  led  to  this  configuration. 

With  any  given  temperature  distribution  in  the  test  section  of  the  apparatus,  that  is, 
the  gradient  along  the  composite  bar  made  up  of  the  specimens  and  the  temperature  variation 
down  the  tantalum  cylinder  supporting  the  guard  heater,  any  lateral  beat  flow  to  or  froqa  the 
specimens  is  proportional  to  the  thermal  conductivity  of  the  Insulation,  h |  .  If  kt  la  an 
average  conductivity  of  the  specimen  composite  bar,  then  the  importance  of  such  spurious 
heat  flows  in  apparatus  performance  diminishes  as  the  ratio  kf  /k(  decreases. 

For  the  alumina  bubble  insulation-beryllium  oxide  combination,  studied  here,  this 
ratio  varied  from  about  0.  0008  at  room  temperature  to  about  0.  02  at  200CT  F.  Simple  heat 
transfer  analysis  of  the  system  showed  that  heat  losses  were  very  small  (on  the  order  of  1 
to  2%  at  most)  even  for  cases  where  there  was  very  poor  matching  between  the  temperature 
gradients  down  the  specimen  cylinder  and  the  tantalum  guard  cylinder. 

This  situation  was  verified  experimentally  as  discussed  in  some  detail  below  In 
connection  with  the  beryllium  oxide  and  tantalum  thermal  conductivity  measurements,  and 
was  also  verified  analytically  with  an  exact  Fourier  analysis  of  the  system  as  discussed 
under  "Heat  Losses.  "  Thus,  the  use  of  the  guard  heater  was  mainly  to  provide  a  nearly 
isothermal  zone  for  specimen  thermocouple  leads  extending  out  into  the  insulation  and  to 
prevent  excessive  power  dissipation  from  the  main  beater,  and  only  secondarily  to  control 
radial  heat  losses  from  the  specimen  column.  Of  course,  if  the  ratio  kj  /k  %  were 
significantly  larger  (e.g. ,  near  unity)  the  above  comments  would  not  be  valid. 

2.  Instrumentation 

In  most  studies  where  conductivity  measurements  are  run  with  cylindrical  specimens, 
the  thermocouples  are  mounted  on  the  specimen  surface.  However,  in  the  case  of  thermal 
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contact  resistance  tests,  it  was  found  from  the  literature  that  most  investigators  mounted  the 
couples  in  radial  holes  drilled  to  the  specimen  centerline.  In  the  brief  analysis  which  follows 
the  variables  determining  the  best  radial  location  of  the  thermocouples  are  discussed  for  the 
particular  system  used  here. 

The  first  limiting  case  considered  is  where  the  thermocouples  are  located  in  radial  holes, 
the  thermal  conductivity  of  the  thermocouple  sensor  kT  ^  being  much  less  than  that  of  the 
specimen.  ks.  This  is  represented  by  the  situation  where  a  couple  in  an  insulated  ceramic 
tube  is  used  in  measurements  with  metals  or  high  conductivity  nonmetalllcs.  The  extent  to 
which  the  couple  perturbs  the  heat  flow  in  the  specimen  can  be  readily  established  using 
equations  analogous  to  those  developed  for  the  fluid  dynamics  problem  of  incompressible, 
irrotatlonal  flow. 


If  a  stream  function  ^  is  defined,  then  the  lines  of  constant  heat  flux  are  represented 
geometrically  by  families  of  curves  where  ^  ■  constant.  In  the  case  of  steady-state 
two-dimensional  heat  transfer  the  Laplace  equation  in  cylindrical  coordinates  can  be  written  as 


J_  d  ^  I  dz\p 

dr*  r  dr  r*  dQ * 


(I) 


The  colutlon  of  this  equation  for  heat  flow  around  a  nonconducting  cylinder  of  radius  a  in  an 
isotropic  media  is 


^  =  <1. 


sin  Qm[2r~  r  ) 


(2) 


The  Important  boundary  condition  which  permits  use  of  this  analytical  approach  for  the  heat 
transfer  case  is  the  condition  that  the  ^  ■  0  streamline  conforms  to  the  surface  of  the 
perturbing  cylinder.  A  plot  of  thin  function  is  given  in  Figure  6. 

The  extent  of  the  perturbation  in  the  direction  Q  *  90*  is  given  by  the  rate  of  change 
of  ^  with  respect  to  r.  The  perturbation  is  zero  when  this  quantity,  |di///dr  j  0  , 

becomes  constant.  From  Equation  2 

s  q  o*  (I  -f  \)  (3) 

ao  r* 

0- 90* 
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Figure  6 
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Thus,  the  perturbation  diminishes  rapidly  with  radial  distance,  and  is  approximately  10%  at 
a  distance  of  r  »  3a  (one  diameter  away  from  the  surface  r  -  a).  As  the  conductivity  of  the 
specimen  approaches  that  of  the  cylinder,  the  perturbation  at  a  given  radial  distance  decreases 
The  approximate  variation  as  a  function  of  specimen  conductivity  is  indicated  as  the  ’'Internal" 
curve  in  Figure  7. 


Now  considering  the  alternate  case  of  surface  mounting  of  the  thermocouples,  major 
errors  could  be  encountered  in  the  temperature  measurement  is  radial  beat  losses  become 
significant.  These  losses  are  proportional  to  kj  (  dj  /dr  )•  If  the  overall  temperature  drop 
down  the  composite  bar  made  up  of  meters  and  specimens  is  assumed  to  be  constant  ihen  the 
radial  heat  loss  will  be  approximately  constant  for  a  given  insulation.  On  the  other  hand,  if 
this  same  overall  drop  is  maintained  from  test  to  test,  the  specimen  heat  flux  will  be  pro¬ 
portional  to  the  mean  conductivity  of  the  composite  bar,  V8.  The  percentage  heat  loss,  which 
is  a  direct  indication  of  possible  specimen  temperature  measurement  errors  for  surface 
mounted  thermocouples,  1b  then  given  by  the  following  equation 


qio»« 

‘’specimen 


D(r, 


z) 


AL 

dr 

dJ 

dz 


(4) 


Thus,  for  a  given  system  the  temperature  measurement  error  is  approximately  proportional 
to  (1  As).  In  Equation  4  conditions  typical  of  those  used  in  this  investigation  were  used  to 
generate  the  "Surface"  curve  shown  in  Figure  7. 


It  is  clearly  shown  that  the  optimum  location  of  the  thermocouple  sensor  for  minimum 
error  in  the  temperature  measurement  varies  as  a  function  of  the  conductivity  of  the  specimen. 
Because  the  thermal  conductivity  of  the  beryllium  oxide  decreases  rapidly  with  temperature 
it  was  concluded,  based  on  the  curves  shown  in  Figure  7,  that  both  surface  and  internal 
thermocouples  should  be  used  in  studying  the  performance  of  this  apparatus.  For  the 
partlrular  experimental  arrangement  used  in  this  investigation  the  curves  are  quantitative. 

For  other  systems  the  curves  Indicate  qualitatively  that  as  the  thermal  conductivity  of  tne 
specimens  under  test  become  lower,  it  may  be  more  desirable  to  mount  the  measurement 
thermocouple  internally. 
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Figure  7.  Temperature  Field  Perturbation  as  a  Function  of  Specimen  Conductivity 


3.  Heat  Loss  Analysis 

In  order  to  arrive  at  an  analytically  exact  prediction  of  the  beat  exchange  between 
the  specimen  column  and  the  insulation  surrounding  it,  it  is  necessary  to  lave  an  expression 
for  the  radial  temperature  gradient  variation  down  the  column  surface.  If  this  is  available 
then  an  equation  of  the  following  form  can  be  applied, 


«lo..  ‘l 


(#> 


The  expression  for  ( dT/dr  )  was  obtained  by  developing  a  general  equation  for  the 
temperature  field  in  the  annular  Insulation  region  surrounding  the  specinlens  and  performing 
the  indicated  differentiation.  This  temperature  field  equation  was  constructed  by  super¬ 
position  of  solutions  for  the  steady-state  boundary  value  system  made  up  of  the  insulation 
region  with  its  sun  ace  temperature  distributions.  Hie  most  general  solution  developed  for 
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thin  axisymmetric  system  was  for  arbitrary  temperature  distributions  fQ  (z),  fj  (z),  f2  (r) 
along  the  inner  and  outer  cylindrical  boundaries  and  the  upper  surface  boundary  respectively. 
The  lower  surface  boundary  was  assumed  to  be  at  a  constant  temperature,  T  =  O,  although 
this  is  essentially  a  temperature  field  normalization  condition  which  could  very  easily  be 
relaxed.  A  sketch  of  the  system  and  these  boundary  conditions  is  given  in  Figure  8. 

A  complete  derivation  of  the  insulation  temperature  field  equation  TlnflUpf(r,  z)  as  well 

as  the  equation  for  qjos8  is  included  in  Appendix  I.  The  discussion  below  summarizes  the 
Important  aspects  of  the  derivation. 

Because  the  governing  differential  equation,  V2  T  =  O,  is  linear  (i.e.,  kj  * 
f(T,  r,  z)),the  various  boundary  conditions  where  considered  separately,  the  solution  for  each 
then  added  to  give  the  required  result.  A  solution  was  developed  for  boundary  conditions  1 
and  2  simultaneously,  and  then  a  solution  for  boundary  condition  4. 


BOUNDARY  CONDITIONS 


(1 ) 

«■='<> 

T  =  f0{  i ) 

(2) 

r  =  r. 

*  , 

T  =  f,{z) 

(3) 

ro^  r  ^r«. 

*  =  o, 

H 

ii 

O 

(4) 

ro  *  r  ^  r,  • 

*  =  L, 

T=  f,(r  ) 

Figure  8.  Insulation  Segment  and  Boundary  Conditions 
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The  solution  of  Laplace's  equation  in  cylindrical  coordinates  for  the  axisymmetric  case 
is  performed  by  separation  of  variables  with  integration  of  the  resulting  total  differential 
equations.  For  the  case  when  the  boundary  conditions  (B»  C. )  are  a  function  of  the  axial 
coordinate  (B.  C.  1  and  2)  the  general  product  solution  is  in  terms  of  ein  {*.)  and  modified 
zero  order  Bessel  functions  of  the  first  and  second  kind,  IQ  (r)  and  KQ  (r).  When  the  boundary 
conditions  are  a  function  of  the  radial  coordinate  (B.  C.  4)  the  general  product  solution  is  in 
terms  of  the  hyperbolic  function  sinh  (z)  and  zero  order  Bessel  functions  of  the  first  and 
second  kind,  JQ  (r)  and  ^  (r). 

Considering  first  the  analysis  involving  boundary  conditions  1  and  2,  separation  of 
variables  leads  to  the  following  general  solution  of  Laplace's  equation* 


T  ( r,*l  *  [a*I0  (  Xr)  +  b'ko  l\r)  j  £  c'cm  (X  z)  +•  0*  sin  (  X*  ) ]  (  6  ) 

where, 

A',  B',  C\  D'  *  constants 

The  T  s  O  boundary  condition  along  the  surface  z  =  O  leads  to  the  requirement  that 

C'  *  0 


The  T  -  O  boundary  condition  along  the  surface  z  =  L  lead  to  the  following  eigenvalue 
restrictions  on  the  parameter  X  . 


Xl_=  n  7T  ,  n  sO,  1 , 2  ,  •  •  • 

Thus,  a  particular  solution  of  the  differential  equation  becomes, 


(7) 


The  constants  in  this  equation  are  now  subscripted  as  E'  and  F'  since  there  is  a  solution 

n  n 

for  each  value  of  n. 


To  satisfy  boundary  conditions  1  and  2  it  is  necessary  to  utilize  an  infinite  seriea 
solution  of  the  differential  equation  constructed  as  a  linear  combination  of  solutions  of  the 
form  given  in  Equation  7.  Then 


T  (r, z  ) 


I  T(r,  z)„ 
n  =  I 


(8) 
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This  linear  combination  is  also  a  solution  of  the  original  differential  equation.  After  a 
regrouping  of  the  constants  in  cylindrical  function  form  the  infinite  series  solution  can  be 
written 


T(r,i)=  I 


^(~Q~Ti~)Gn(r  »f o)  +  ^>(iLT&~)Gn(r.»  r  [ 
6  (r,;r2) 


(r\m\ 
i,n  (—) 


(9) 


where 


(H  7T  Tq  y  /  /  n  7T  Tq  \ 

-r2") ■  E"  xo  (— r2-) *  F, 


n  7T  r0 


>/  /  n  7T  fp 

«  v  L 


) 


nirr,  ■ 

\  /_ n7r  r0  > 

\  -T  1 

f  n  7T  rQ  > 

\  K  (n1rf\ 

L 

/  Ko  \  L  t 

I  Io  ( 

i  L  1 

i  Ko(— r 

Substituting  B.  C.  1  and  2  in  the  above  equation  and  developing  the  Fourier  coefficient 
integral  conditions  in  the  usual  way  for  this  orthogonal  family  (Reference  2)  leads  to  the 
following  general  solution  of  the  differential  equation  meeting  boundary  conditions  1  and  2. 


n /  n7rz  \  L 

=  0nCr,;rol  {G",,;,o1/  dz 

+  Gn  (r,‘  r  >  /  f0  (z)  5,n  (— j7")  dz  } 


(10) 


Considering  now  the  solution  involving  boundary  condition  4,  separation  of  variables 
leads  to  the  following  general  solution  of  Laplace's  equation: 


T(r,i  )=  [a  J0  l^r)+BY0Kr)  ]  [c  cosh  (  /z  )  +  D  tinh  (ID 

where 

A,  B,  C,  D  ■  constants 

The  T  =  O  boundary  condition  along  the  surface  z  =  O  leads  to  the  requirement  that 
C  «  O 

Applying  the  boundary  conditions 


r  =  r0  ,  z  .  T  :  0 

r  :  .  1  ,  T  :  0 
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and  eliminating  the  constants  A  and  B  from  the  resultant  equations  leads  to  the  following 
condition  which  must  be  fulfilled: 

Yo  ( K  ro  >  Jo  f .  )  ~  Yo  < <  >  J0  ( K  ro  >  *  0  < 

The  above  condition  is  met  by  an  infinite  number  of  eigenvalues,  7^,  which  are  the 
positive  roots  of  this  equation.  Introducing  cylindrical  function  notation, 

uo  (  K n  ’  ro  •  ri  >  *  Jo  (  Kn  ro  )  Yo  '•  )  ~  Jo  (  Vnri  1  Yo  {Kn' o  > 

We  can  now  define  the  eigenvalues,  ^n,  by 

uo  C<„;  ro  .  r,  >  =  O  ,  n  =  1 , 2,  3  , .  I 

Thus,  a  particular  solution  of  the  differential  equation  becomes 


T(».«>  =  [*„  J0  <*■„'  I  +  B„  Y0  <*„o][d„  (15) 

where  again  the  constants  are  subscripted  to  reflect  the  family  of  solutions  represented  by 
the  eigenvalues,  n. 

Writing  the  corresponding  infinite  series  solution  and  recasting  the  form  of  the  constants 
leads  to  the  following  representation  for  T(r,  z), 


>o«(n  > 


wliere 


<2  <  '  >*  En  Jo  lXnr  '  +  F„\  <*„'> 


Substituting  B.  C.  4.  in  the  above  equation  and  developing  the  Fourier  coefficient  integral 
conditions  leads  to  the  general  solution  given  below  in  Equation  17  as  a  cylindrical  function 
Fourier  series.  The  norm  N  |uQ(r)  }  for  this  class  of  orthogonal  functions  is  found  by 
applying  Lommel  integral  relations. 


T(  ,  -  77,2  y  Kn  Jo  *Kn  ro*  uo^Kn  •  r  •  ri  *  sinh  Knz  f  1  .  .  .  .  .  . 

Tlr.zl-  -yI  — - - - : - J  7  ft<r)«0<Kn*r*r.,dr  (,7> 

0  =  1  [J0<Kn'o>  -Jo<Knr,  )J*inh  KnL  fo 
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Hie  general  solution  of  the  original  differential  equation  meeting  all  tiv*  boundary 
conditions  given  in  Figure  8  is  simply  the  sum  of  Equations  10  and  17. 

For  this  particular  experimental  arrangement,  boundary  condition  4  is  not  arbitrary. 

The  constraint  oh  f2(r)  is  of  a  physical  nature  rather  than  mathematical.  If  fj(L)  and  fQ(L) 

are  the  temperatures  on  the  z  *  L  plane  at  points  r  =  r^  and  r  -  r  ,  respectively,  then, 
as  developed  in  Appendix  I,  to  satisfy  y  2T  -  O, 

JL 

T.I0(L)+{|,  (L)-(0(L)}[ - '-f—  ]  118) 

It  will  be  shown  later  than  fQ<z)  and  fj(z)  can  be  quite  accurately  represented  in  the  form, 

f0  (s  )  =  T0  ( I -a )  +  a  TLo  (19) 


where 


f ,  (r  )  =  T,  (  Ha  )+  a  Tl 


(20) 


• =  ir 

'o  =  'oto).rLo»  f0(L) 


Using  Equations  18  -  20  for  explicit  evaluation  of  the  various  Fourier  coefficient  integrals 
leads  to  the  following  equation  for  the  temperature  distribution  in  the  insulation  segment- 


T(r. 


a 

m 

.  y 

r  ^  Kn  fo  ^uo  •  r  -ro 

)sinh  (Knr  )  j 

M 

^  V 

z. 

|  n=l 

L 

[•>1  -j'o<Xn  'll 

sin h  (Kn  L  )  J 

lk„  (l+^r0)-T  l'o>}  {j0‘Kn'o»- Jo  <«nf,  »} 


ITU|-TLolk'K»'ol-^ri  ~^o fo } ] 


+  L  C"|7v|{G"lf'r»l  +  •■'l(^)[VTLol-o"]}  12 


I) 


22 


AFML-TR-69-1 


Equation  5  is  now  applied  to  determine  the  total  heat  exchange  between  the  composite 

specimen  bar  and  the  insulation.  In  using  this  equation  we  assumed  that  the  insulation  thermal 

conductivity,  Kj,  was  temperature  and  position  independent.  A  temperature  and  spatially 

dependent  conductivity  could  be  introduced  in  evaluating  the  Equation  5  integral  although  it  is 

important  to  bear  in  mind  that  this  is  still  only  an  approximation  since  the  use  of  the  condition 
2 

V  T  *  O  implies  temperature  independent  properties. 


The  final  result  of  the  substitution  and  integration  of  Equation  5  is  as  follows. 


-  2  rr  r.  k 

f(dM  •• 

loss  o  1  J 

( 

3  v  ar  jr  =  ro 

27T.0k|  {2  I 
n=l 

rGn(r9'r<?)  /T  T  (|)n\t  6«(ro.ri>  /, 

L  Gn(r,  ;rQ)  lT«  V  °  r  Gn(r,,r0)tTo 

2 

"  r  _  Jo(Knr.)  r  COSh(KnL)-l 

'Mr) 

... €"  *  *<w  L 

]} 


(22) 


where 


f„  =  VKnro|/[{Jo<Kn'o>-Jo,(lV.>}  •i"h  <Kr,L>] 

Vn  =  tl0  *  [Jo  ^  Knro*  —  ’VKi>ri  *  ] 

+  (T  _Ti  >  {Jo<Knro>  A*  ri  -  J0<  Knri  )  A*  'o  } 

L,  V  I 


/  .  /n7rr\  / n  irrn  \  /  n  7T  rn  \  /  n  trr  \ 

°„,r -'o'- 1. (— )*„(-rL)+  ^(-T-2-)  *,(—[—) 


Alternate  terms  of  the  first  summation  are  zero  due  to  the  factor  |— (-  |  )nj;  hence 
only  odd  values  of  n  need  be  considered.  In  the  second  summation,  a  determination  of  the 
eigenvalues  as  defined  by  Equation  14  is  necessary.  An  infinite  series 

approximation  given  by  Gray  and  Mathews  (Reference  3)  was  applied;  limited  tabulations  are 
also  found  in  Jahnke  and  Emde  (Reference  4).  The  first  few  roots  are  given  in  Table  I  for  the 
case  where  r  j/rQ  =  6.0.  For  n  >  6  the  roots  are  given  to  within  0. 1%  by  Equation  23, 


(23) 
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TABLE  I 

POSITIVE  ROOTS  OF  u^jf'  ;  rQ;  *  0  WHEN  rj/r^  *6.0 


n 

ro  fn 

n 

ro  ’tn 

;  i 

0.  226 

5 

3. 134 

2 

1.233 

6 

3.764 

3 

1.874 

7 

4.393 

4 

2.  506 

8 

5.023 

Digital  computer  evaluation  of  Equation  22,  a  discussion  of  its  convergence  characteristics, 
and  a  comparison  of  the  computer  values  with  experimental  heat  loss  measurements  are 
given  in  the  following  section,  which  concerns  the  experimental  performance  of  the  apparatus. 


C.  APPARATUS  PERFORMANCE 

1.  Thermal  Conductivity  Specimens 

For  the  thermal  contact  resistance  measurements,  it  is  necessary  to  know  the  thermal 
conductivity  of  the  test  specimens  quite  accurately  since  these  specimens  are  used  as  heat 
meters.  A  variety  of  thermal  conductivity  measurements  were  conducted  on  beryllium  oxide 
(BeO)  and  tantalum  with  two  objectives  in  mind:  (1)  to  assess  apparatus  performance,  and 
(2)  to  obtain  accurate  thermal  conductivity  data  on  the  BeO  for  the  later  use  in  the  contact 
resistance  measurement  series. 

The  beryllium  oxide  specimens  were  obtained  from  Coors  Porcelain  Company, 

Golden,  Colorado.  The  Coors  body  designation  was  BD-995;  the  composition  obtained  by 
spectrographic  analysis  is  given  in  Table  II.  The  specimens  were  isostatically  pressed  from 
powder  stock  with  average  particle  size  of  approximately  one  micron;  prefired  to  1100*  C 
followed  by  drilling  of  the  thermocouple  instrumentation  holes;  and  finally  fired  at  1690*  C  for 
approximately  4  hours.  Based  on  a  theoretical  density  of  3.01  the  specimens  were  94.7%  dense. 
On  the  conductivity  specimens  the  flat  end  faces  were  lapped  to  within  3  fringes  (about  35  /xin. ). 

The  tantalum  specimens  were  obtained  from  Fansteel  Metallurgical  Corp. , 

Chicago,  Illinois.  The  sintered  rod  stock  was  unannealed,  metallurgical  grade  and  was 
fabricated  into  specimens  in  this  condition.  The  composition  is  given  in  Table  III.  The  theoretical 
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TABLE  n 

COMPOSITION  OF  BERYLLIUM  OXIDE  TEST  SPECIMENS 


CONSTITUENT 

WEIGHT  PERCENT 

A12  °3 

0.30 

sio2 

0.18 

MgO 

0.10 

NazO 

0.05 

CaO 

0.  02 

Fe2°3 

0.01 

BeO 

99.34 

3 

Density:  2.  85  gms/cm 

Geometry:  Cylinders  -  1. 0  inch  diameter  x  1.  5  inch  length 
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TABLE  III 

COMPOSITION  OF  TANTALUM  TEST  SPECIMENS 


CONSTITUENT 

WEIGHT  PERCENT 

Ta 

99.  8  Minimum 

Cb 

0.02 

C 

0.03 

Ni 

0.  02 

Fe 

0.  02 

Ti 

0.01 

Si 

0.  02 

W 

0.  03 

H 

0.01 

N 

0.015 

O 

0.  03 

Density: 

Geometry: 

3 

16.  56  gms/cm 

Cylinders  -1.0  inch  diameter  x  2.  0  inch  length 

3 

density  was  near  16.66  gms/cm  ,  thus  this  specimen  material  was  99.4%  of  theoretical.  As 

a  check  of  the  metallurgical  condition  of  the  specimens ,  extensive  Vickers  hardness 

measurements  were  made.  Generally,  the  Vickers  hardness  of  tantalum  is  approximately 

160  for  unannealed  material  and  drops  to  around  90  for  fully  annealed  metal.  Measured 

2 

values  of  our  tantalum  specimens,  In  units  of  kgms/mm  ,  were  nearly  constant  at  125, 
Independent  of  load  level  variations  from  150  to  1000  grams. 
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2.  Thermal  Conductivity  Measurements  on  BeO 

For  the  beryllium  oxide  thermal  conductivity  determinations,  a  1.0  inch  diameter 
x  1.  5  inch  length  specimen  was  positioned  between  two  Armo  Iron  cylinders  of  the  same 
diameter.  Armco  iron  is  a  widely  used  reference  standard  in  comparative  measurements, 
its  conductivity  being  known  with  uigh  accuracy  (Reference  5).  Very  thin,  0.  0001  inch  platinum 
foil  was  placed  between  the  contacting  surfaces  to  reduce  thermal  contact  resistance. 

Reference  grade  0.008  inch  diameter  Pt/Pt-10  Rh  thermocouples  were  used 
throughout.  The  Armco  Iron  specimens  had  couples  mounted  in  grooves  on  the  surface  and 
couples  inserted  radially  to  the  specimen  centerline.  The  centerline  couples  were  mounted 
in  two-bore  42  mil  OD  alumina  tubing  (high  purity  stock  from  McDanel).  In  all  cases  the 
thermocouple  leads  extended  radially  in  a  nearly  isothermal  plane  for  a  distance  of  at  least 
2  inches.  Temperature  readings  were  taken  with  a  Leeds  and  Northrup  K-3  potentiometer. 
Thermocouple  calibrations  were  run  before  and  after  the  conductivity  determinations. 

Maximum  deviations  were  0.2%.  All  measurements  were  conducted  in  vacuum,  the  level 
varying  between  1.  0xl0-5  and  1.  2xl0-6  torr. 

Two  typical  temperature  distribution  plots  during  the  conductivity  determinations 
are  given  in  Figure  9.  A  dead  weight  compressive  load  of  375  psi  was  maintained  on  the 
test  column  to  reduce  interfacial  temperature  drops.  The  average  contact  conductance  of  the 
BeO/Armco  interfaces  with  the  platinum  foil  shim  was  about  1000  BTU/hr-ft2  -  *  F.  At  a 
375  psi  load  level  this  is  a  fairly  high  value  for  an  interface  with  a  brittle  material  component, 
and  indicated  that  the  platinum  foil  was  effective.  Figure  9  shows  the  consistency  of  the 
temperature  readings  to  be  quite  good.  The  guard  was  adjusted  such  that  it  matched  closely 
the  temperature  gradient  in  the  BeO  test  specimen. 

The  thermal  conductivity  results  are  given  in  Figure  10  along  with  several  curves 
from  the  literature  for  BeO  specimens  in  the  same  density  and  purity  range  (References  6-9). 
For  clarity,  a  curve  has  not  been  drawn  through  the  data  points  obtained  in  this  investigation. 
The  precision  of  the  data  was  found  to  be  better  than  ±2%.  An  exact  accuracy  figure  relative 
to  literature  results  was  not  calculated  since  slightly  differing  curves  were  obtained  by 
\arious  investigators  using  material  of  the  same  composition  studied  here.  A  3%  scatter 
interval  is  shown,  hovever,  to  indicate  the  close  agreement  between  literature  values  and 
the  data  obtained  here.  The  experimental  data  are  tabulated  in  Appendix  IV. 
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Figure  9.  Typical  Temperature  Distributions  During  BeO  Thermal  Conductivity 
Measurements 
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Figure  10.  Thermal  Conductivity  of  Beryllium  Oxide  as  a  Function  of  Temperature 
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For  beryllium  oxide  and  other  dielectric  materials  the  thermal  conduction  process  is  due 
entirely  to  phonon  transport,  the  electronic  contribution  being  negligible.  Applying  the 
conventional  mean  free  path  concept,  this  lattice  thermal  conductivity  can  be  expressed  as 

ki  :  T  Cv  V<£  (24  ) 

Correlations  of  phonon  scattering  predict  that  for  crystalline  dielectrics  is  nearly 
proportional  to  1/T.  Since  V  is  a  weak  function  of  temperature  (proportional  to  the  0.  5  power 
of  the  elastic  modulus)  and  since  Cv  is  approximately  constant  in  the  temperature  range 
studied  here,  it  is  expected  that  the  thermal  conductivity  of  BeO  would  be  proportional  to  1/T. 
This  dependence  should  be  obeyed  to  temperatures  of  about  3000*  F,  where  radiation  heat 
transfer  starts  to  become  significant  (Reference  10). 

The  experimental  results  of  this  investigation  are  replotted  in  this  form  in  Figure  11 
together  with  data  from  the  literature  (References  6,  7,  and  11).  The  "recommended"  curve 
in  this  figure  was  formulated  by  the  Thermonhysical  Property  Research  Center  and  based 
on  an  analysis  of  about  50  different  sets  of  literature  data  on  BeO.  This  curve  is  for  material 
of  theoretical  density,  P  ■  3.01.  Its  agreement  wi‘h  theory  is  excellent,  the  data  clearly 
defining  a  straight  line  Further,  this  line  falls  between  the  curves  obtained  by  other 
invesUii-itors  for  material  in  the  same  density  and  purity  range.  The  experimental  data 
pole's  Powell  and  Day  are  not  dotted  lr  Figure  11  to  avoid  cluttering.  Separate  plots  of 
thei.’  daU  produced  scatter  about  the  same  as  that  observed  in  this  investigation.  Based  on 
this  correlation  the  conductivity  data  from  this  investigation  wus  extrapolated  to  1800*  F,  this 
being  the  maximum  temperature  to  which  conductivity  data  were  required  in  the  contact 
resistance  computations.  The  highest  temperature  data  point  wtis  taken  at  1120*  F. 

3.  Thermal  Conductivity  Measurements  on  Tantalum 

Measurements  with  tantalum  were  performed  in  a  manner  identical  to  those  with  BeO. 
The  scatter  of  the  data  points  about  a  mean  line  through  them  was  a  small  degree  higher  than 
that  observed  with  the  BeO,  the  maximum  amounting  to  *1.  59L  The  data  are  plotted  in 
Figure  12  together  with  data  from  the  literature.  The  numerical  value  fabulation  is  included 
in  Appendix  D. 

The  most  probable  value  of  thermal  conductivity  for  this  metal  is  more  open  to  question 
than  for  the  beryllium  oxide,  as  is  apparent  from  the  curves  in  Figure  12.  For  this  reason 
the  recommended  curve  is  somewhat  uncertain.  The  data  from  this  work  agree  with  the 
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Figure  12.  Thermal  Conductivity  of  Tantalum  as  a  Function  of  Temperature 
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recommended  curve  in  the  range  from  room  temperature  to  50(f  F,  but  because  of  a  difference 
in  slope,  there  is  a  7.  5%  disagreement  at  the  highest  test  temperature,  115(f  F.  However, 
other  literature  deta  at  much  higher  temperatures  show  disagreement  with  the  recommended 
curve  ranging  from  0. 1  to  15%,  and  the  deviations  are  in  the  same  direction  as  those  for  the 
data  obtained  here.  From  limited  electrical  resistivity  measurements  and  application  of  the 
Wiedeman- Franz  relationship  it  has  been  estimated  that  the  electronic  component  of  the 
thermal  conductivity  of  tantalum  is  approximately  88%  (Reference  12).  The  lattice  component 
is  thus  about  12%.  Some  oi  the  disagreement  among  authors  indicated  in  Figure  12  could  be 
resolved  if  electrical  resistivity  data  were  available.  Such  measurements  have  not  yet  been 
made  on  the  tantalum  studied  in  this  investigation  since  it  was  beyond  the  scope  of  the  project. 
It  is  felt,  however,  that  the  tantalum  conductivity  results  are  a  further  demonstration  of  the 
successful  performance  of  the  equipment  used  in  this  report. 

4.  Heat  Loss  Measurements 

During  the  beryllium  oxide  conductivity  determinations,  data  waB  obtained  on  heat 
exchange  between  the  insulation  and  the  specimen  column.  This  data  consisted  of  measured 
differences  between  the  heat  flux  through  each  of  the  meter  bars  as  a  function  of  absolute 
temperature  level  and  variations  of  the  flux  in  the  specimen  column  as  the  guard  heater 
temperature  distribution  was  changed. 

To  evaluate  Equation  22,  it  was  necessary  to  have  vacuum  environment  values  of  the 
alumina  bubble  insulation  conductivity,  kj,  as  a  function  of  temperature.  The  rationale  for 
determining  this  functionally  is  given  in  Appendix  III.  Table  11,  Appendix  A,  gives  a  summary 
of  the  temperature  distribution  and  insulation  conductivity  data  required  in  Equation  22. 

In  evaluating  the  series  in  Equation  22,  term  by  term  printout  of  the  computer  results 
was  employed  to  observe  convergence  characteristics.  Hand  calculation  of  several  terms 
was  also  run  to  check  the  computer  results.  The  first  series  was  carried  to  N  -  40. 

For  n  >  20  each  term  contributed  no  more  than  0.3%  to  the  sum.  In  the  case  of  the  second 
series,  the  first  term  in  the  summation,  represented  by  Kj,  varied  between  70  and  93%  of 
the  series  sum.  Subsequent  terms  alternated  sign  with  decreasing  magnitude  as  might  be 
expected  for  a  series  having  linear  JQ  (Knr)  factors.  Summation  through  was  performed 
to  obtain  a  sum  varying  less  than  1%  with  inclusion  of  higher  terms. 

The  comparisons  between  the  experimental  heat  loss  measurements  and  those 
predicted  bv  Equation  22  are  given  in  Table  IV  as  percentages  of  the  average  heat  flux  in  the 
specimen  column.  Positive  values  represent  heat  flow  from  the  specimen  column  into  the 
insulation. 
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TABLE  IV 

EXPERIMENTAL-ANALYTICAL  HEAT  EXCHANGE  COMPARISONS 


Data  Point 

Average  Heat  Flux  in 
the  specimen  Column 
(BTU/hr) 

Predicted 

Value 

% 

Experimental 

Value 

% 

2 

37.4 

2.6 

5.4 

4 

62.7 

2.8 

6.4 

9 

193.  5 

7.  1 

5.  5 

10 

230.7 

9.  0 

5.4 

11 

230.3 

10.0 

5.4 

12 

229.3 

11.9 

4.  1 

13 

276.  5 

14.0 

4.0 

15 

354  .0 

16.8 

-1.  1 

18 

443.7 

19.7 

-5.  5 

If  is  tin?  midpoint  temperature  of  the  beryllium  oxide  specimen  and  T  ar^ 

is  the  guard  heater  temperature  at  the  same  axial  position,  then  (YgUar(j/(TgeQ)m  -  1)100 
is  the  percentage  by  which  the  guard  temperature  exceeds  that  of  theBeO  specimen;  this 
quantity  will  be  termed  the  guard  mismatch.  For  data  points  10-12,  the  mean  temperature  of 
the  BeO  was  held  approximately  constant,  while  the  guard  mismatch  varied  from  8  to  23%. 

A  .een  in  Table  IV  the  heat  exchange  between  specimen  column  and  insulation  varied  only 
1.3%  as  observed  experimentally  and  only  2.9%  according  to  the  analytical  prediction.  Thus, 
thermal  guarding  of  the  specimen  column  does  have  an  effect  on  the  heat  flowing  in  the  column 
but  as  suggested  in  the  design  philosophy,  the  effect  is  not  of  high  sensitivity.  The 
experimental  conductivity  data  was  obtained  with  guard  mismatching  ranging  from  +23  to  -11%. 
Across  this  wide  range  the  average  heat  exchange  between  insulation  and  specimen  column 
was  4. 7%. 

It  is  important  to  note  that  there  is  some  uncertainty  in  the  experimental  heat  exchange 
values.  For  data  points  10  and  above  these  exchange  values  are  obtained  bv  taking  the 
difference  between  large  heat  flux  values,  hence  uncertainty  in  these  fluxes  greatly 
accentuates  uncertainties  in  the  heat  exchange  values.  For  data  points  1  to  9  the  principal 
uncertainty  in  the  heat  exchange  values  is  the  error  in  the  temperature  gradient  determinations 
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in  the  meter  bars  above  and  below  the  BeO  specimen.  This  uncertainty  is  particularly 
important  for  the  lower  temperature  data  points  where  the  overall  temperature  drops  in  the 
meter  bars  are  relatively  low. 

The  change  in  sign  of  the  experimental  heat  exchange  values  for  the  higher  temperature 
data  points  are  internally  consistent.  That  is,  in  the  cases  where  the  exchange  is  positive 
(heat  flow  into  the  insulation),  the  thermocouples  on  the  specimen  surface  give  lower  readings 
than  those  in  the  same  plane  at  the  specimen  centerline.  When  the  experimentally  observed 
heat  exchange  values  change  sign,  the  surface  thermocouples  are  found  to  read  higher  than 
the  corresponding  centerline  couples. 

Now  a  comparison  of  the  experimental  heat  exchange  values  with  the  analytical 
predictions,  based  on  the  uncertainties  in  the  experimental  values  discussed  above,  shows 
the  predictions  consistent  at  least  through  data  point  13. 

The  last  two  data  point  comparisons  in  Table  IV  require  further  elaboration.  The 
experimental  thermocouple  data  reveals  an  uncertainty  of  gradient  in  the  meter  bars  of 
3.  1%.  In  the  cases  where  the  measured  heat  exchange  is  near  zero  and  going  negative  this 
uncertainty  in  the  meter  bar  gradient  introduces  a  particularly  large  uncertainty  in  the  heat 
exchange  value.  The  range  for  data  point  I Is  from  -1.1  to  +2.8%,  and  for  data  point  18, 
from  -5.  5  to  -2.  2%. 


Further,  the  magnitude  of  the  h(  at  exchange  is  large  here  compared  with  the  low 
temperature  cases.  For  data  points  15  and  18  this  radial  flux  is  from  10  to  12  times  that 
for  data  point  2.  For  these  higher  temperature  data  points,  the  comparatively  high  radial 
fluxes  introduces  the  possibility  of  significant  temperature  drops  at  the  specimen-insulation 
interfaces  and  the  insulation-guard  heater  interface.  These  temperature  discontinuities  result 
from  the  contact  resistance  associated  with  the  interfaces  and  become  significant  when  the 
radial  flux  is  large.  If  a  conservatively  high  interface  thermal  contact  conductance  of 
20  BTU/hr-ft  -*  F  is  assumed,  the  interfacial  temperature  drop  is  at  least  20*  F  for  the  last 
two  data  points  in  Table  IV.  This  factor  is  significant  with  respect  to  the  analytical 
predictions  because  the  temperature  distributions  used  in  evaluating  Equation  22  were  not 
measured  directly  in  the  insulation  but  in  the  specimen  column  and  the  guard  heater  ring. 

This  is  certainly  one  of  the  largest  factors  in  the  disagreement  between  the  heat  exchange 
values  for  the  high  temperature  data  points. 
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In  general,  the  heat  exchange  predictions  are  of  the  right  magnitude.  However,  it  is 
apparent  that  to  provide  a  more  rigorous  test  of  heat  exchange  predictions  it  would  be  necessary 
to  conduct  extensive  measurements  with  that  specific  aim  in  mind.  Such  measurements  would 
be  most  difficult  to  perform  due  to  the  necessity  of  accurately  measuring  temperatures 
directly  in  the  low  density  insulation  material.  The  temperature  field  equations  developed 
earlier  can  easily  be  generalized  for  predicting  heat  flows  in  any  right-circular  element  with 
arbitrary  boundary  conditions.  For  this  more  general  application,  the  temperature  would  be 
written  in  dimensionless  form  such  as  an  unaccomplished  temperature  ratio,  and  the  physical 
dimensions  of  the  hollow  cylindrical  element  would  also  be  written  in  dimensionless  form. 
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SECTION  m 

CONTACT  RESISTANCE  RESULTS 

A.  CONTACT  SURFACE  CHARACTERIZATION 

Because  contact  surface  condition  has  a  great  influence  on  interfacial  heat  transfer, 
surface  preparation  and  characterization  were  carefully  controlled  from  the  inception  of  the 
program.  It  was  concluded  that  surfaces  with  random  roughness  were  most  desirable  for 
study.  In  this  context  "random"  implies  a  roughness  independent  of  position  on  the  surface. 

To  develop  this  condition,  the  beryllium  oxide  surfaces  were  first  polished  flat  with  diamond 
paste  to  within  3  interference  fringes.  This  was  followed  by  a  vapor-blasting  operation  which 
produced  finishes  in  two  roughness  ranges  nominally  10  and  40 ^  in.  centerline -average  (CLA). 
These  designations  were  used  simply  for  shop  control  of  the  roughness  range. 

Detailed  quantitative  surface  characterization  was  accomplished  with  a  high  precision 
profilometer.  The  instrument  employed  was  a  Taylor-Hobson  "Talysurf"  Model  3  profilometer 
with  a  0.0001  inch  radius  diamond  stylus.  Horizontal  drive  speeds  giving  20  and  100X 
magnification  were  used.  Vertical  amplifications  of  from  1000  to  50, 0O0X  were  used 
depending  on  the  type  of  surface  irregularity  being  studied.  The  physical  configuration  of 
the  varfous  components  is  shown  in  Figure  13. 

Two  types  of  surface  characterization  data  were  obtained:  waviness  and  roughness  traces. 
In  the  waviness  measurements,  relatively  large  scale  surface  undulations  were  recorded  by 
tracing  the  specimen  surface  with  the  stylus  support  arm  moving  along  an  external  optical 
flat.  The  flat  thus  provided  a  reference  line  external  to  the  specimen  surface.  Vertical 
magnification  on  these  traces  was  fairly  low  to  follow  surface  irregularities  of  large  amplitude. 
In  the  roughness  measurements  the  vertical  amplification  was  increased  to  observe  small 
scale  surface  irregularities  having  amplitudes  on  the  order  of  a  few  micro-inches.  To 
prevent  the  surface  waviness  patterns  from  driving  the  roughness  traces  off  scale  at  these 
high  amplifications,  an  external  reference  line  was  not  used.  Rather,  the  specimen  surface 
itself  was  used  as  a  reftrence,  the  stylus  following  a  small  skid  or  shoe  that  moved  ahead 
of  the  recording  tip. 

The  difference  in  surface  character  between  the  10  and  40  p  in.  specimens  is  shown  in 
Figure  14.  Note  that  the  vertical  magnification  on  the  10  ^.in.  specimen  trace  is  twice  that 
of  the  40  jiin.  specimen  trace.  The  10  pin.  trace  reflects  the  higher  degree  of  surface 
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Figure  13.  Surface  Characterization  Equipment 
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polish  used  to  obtain  the  smoother  finish.  For  both  roughness  and  waviness  characterization, 
numerous  traces  were  run  at  different  orientations  across  the  specimen  surfaces.  These 
studies  confirmed  the  random  lay  of  the  surface  roughness. 


A  typical  waviness  trace  is  shown  in  Figure  15.  When  diametrical  traces  of  this  form 
were  run,  it  was  observed  that  the  shape  of  the  profile  was  independent  of  the  orientation  of 
the  specimen.  This  leads  to  the  conclusion  that  the  waviness  irregularities  were  in  the  form 
of  concentric  circular  patterns.  This  observation  is  Consistent  with  the  specimen  fabrication 
and  surface  finishing  procedures  used  and  was  of  considerable  consequence  in  the  formulation 
of  the  interfacial  contact  model  outlined  in  the  discussion  of  results  (Section  IV). 
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Figure  15.  Surface  Waviness  Profile  Trace 


For  the  contact  resistance  correlations  as  a  function  of  compressive  ad,  it  was 
necessary  to  know  the  variation  in  the  sol  id -sol  id  contact  area  as  the  surfaces  were  pressed 
together.  Plots  of  this  form  are  given  in  Figure  16  which  were  generated  in  the  following 
fashion.  First,  a  linear  least  squares  line  was  fitted  through  the  profile  traces  in  Figure  14 
such  that  the  average  distance  between  "peaks"  and  "valleys"  on  either  side  of  the  line  was 
the  same.  A  reference  line  with  the  same  slope  was  then  displaced  upward  until  it  Just 
touched  the  uppermost  peak.  A  family  of  parallel  lines  was  then  drawn  at  various  distances 
between  the  reference  line  and  the  leadt-squares  line.  The  percentage  of  solid  area 
intercepted  by  each  of  these  lines  was  plotted  as  a  function  of  their  distance  from  the 
reference  line  to  give  the  roughness  curves  in  Figure  16.  The  waviness  plots  in  this  figure 
were  generated  in  a  similar  fashion  using  data  of  the  form  shown  in  Figure  15. 
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01  STANCE  FROM  REFERENCE  LINE  -Mi«v 
Figure  16.  Solid  Area  as  a  Function  of  Vertical  Displacement 


/ 


AFML-TR-69-1 

The  quantitative  data  in  Figure  16  reflect  the  surface  characteristics  which  are  evident 
in  Figure  14.  For  the  rougher  surface,  the  increase  in  solid  area  with  vertical  displacement 

is  quite  gradual.  On  the  same  displacement  scale,  the  increase  in  intercepted  area  is  much 
more  abrupt  ror  the  smoother  surfaces.  Even  though  the  waviness  irregularities  are  of 
much  lower  frequency  than  that  of  the  roughness,  the  amplitudes  are  of  the  same  order  of 
magnitude  in  each  surface  finish  range.  Further  discussion  of  these  distinctive 
characteristics  will  be  given  in  Section  rv. 

B.  MEASUREMENT  PROCEDURE 

The  beryllium  oxide  thermal  contact  resistance  specimens  were  of  the  same  configruation 
as  the  thermal  conductivity  specimens:  right  circular  cylinders  of  1. 0  inch  diameter  and 
1.5'inch  length.  Each  specimen  was  instrumented  with  three  centerline  thermocouples 
mounted  at  axial  distances  of  0.20,  0.70,  and  1.20  inches  from  the  contact  interface.  The 
thermocouples,  with  inert  gas  welded  junctions,  were  formed  of  8  mil  pt  versus  8  mil  pt/ 

10%  Rh  and  inserted  in  0.  042  inch  OD  two  bore  tubing  of  high  purity  aluminum  oxide.  To 
ensure  good  physical  contact  with  the  specimens,  the  tips  of  the  couples  were  wound  with  fine 
platinum  wire,  and  the  aluminum  oxide  probes  were  then  force-fit  into  the  0.  048  inch  radial 
holes  drilled  to  the  specimen  centerline.  The  centerline  thermocouple  positions  were 
determined  by  observing  the  inclination  of  rigid  steel  pins  inserted  full  depth  into  the 
instrumentation  holes.  A  traveling  stage  microscope  was  used  for  thiB  purpose. 

After  assembly  of  the  apparatus  in  a  configuration  similar  to  that  shown  in  Figure  2, 
the  thermocouples  were  checked  by  taking  room  temperature  readings  with  respect  to  the 
ice  point.  Post  test  elevated  temperature  calibrations  were  run  also;  in  both  cases  consistent 
readings  within  ±1  volt  (±0.  2*  F)  were  obtained.  A  predetermined  compressive  load  was 
applied  and  after  reaching  a  vacuum  level  of  at  least  5x10  torr,  the  power  to  the  main 
heater  was  adjusted  to  the  desired  mean  interface  temperature .  The  guard  heater  was  then 
adjusted  to  match  the  gradient  in  the  contact  resistance  column.  Typical  temperature 
distributions  were  similar  to  t&ose  shown  in  Figure  9,  except  that  only  beryllium  oxide 
specimens  were  involved.  Depending  on  the  temperature  level,  a  period  of  between  8  and 
24  hours  was  allowed  for  system  stabilization  before  temperature  readings  were  recorded. 

The  longer  timb^periodfl.  were  required  at  low  temperatures  due  to  system  thermal  inertia. 
Provisions  were  made  for  automatic  temperature  adjustment  of  the  guard  heater  relative  to 
the  main  heater;  however,  manual  adjustments  were  used  most  of  the  time  because  the 
thermal  inertia  of  the  system  changed  with  temperature  level.  Steady-state  conditions  were 
defined  such  that,  with  matched  guarding,  the  temperature  variations  were  lep^  than  1.0  ^ 
t  (0.  2*  F)  over  the  period  of  readout  (approximately  1  hour). 
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Temperature  readings  were  taken  with  a  Leads  and  Northrup  K-3  potentiometer  with 
reference  junction  compensation.  Readings  were  made  to  within  ±0. 1  fi  (±0.  02*  F).  The 
magnitude  of  the  interfacial  temperature  drop  was  determined  by  extrapolating  the 
temperature  profiles  in  each  specimen  to  the  contact  interface.  Deviations  of  the  specimen 
gradients  from  linearity  due  to  the  temperature  variation  of  the  thermal  conductivity  had 
an  extremely  small  influence  on  the  calculated  interfacial  drop.  Because  the  overall 
temperature  drops  in  the  specimens  were  small,  no  deviations  from  the  linear  gradient 
condition  were  noted  except  at  the  highest  test  temperature.  To  minimize  errors,  all 
temperature  discontinuity  calculations  were  made  numerically;  graphical  plots  of  the 
temperature  distributions  were  used  only  to  check  relative  consistency  of  the  readings. 

The  thermal  conductivity  data  obtained  in  Section  II.  C.  2.  and  the  measured  temperature 
gradients  in  the  two  BeO  specimens  were  used  to  calculate  the  heat  flux  crossing  the  interface. 
The  ratio  of  the  interfacial  temperature  drop  to  this  average  heat  flux  gave  the  thermal 
contact  resistance  in  units  of  (BTU/hr-ft2-*  F)_1.  The  difference  in  the  heat  flux  in  the 
specimen  above  and  below  the  interface  gave  a  quantitative  measure  of  heat  exchange  with 
the  surrounding  insulation. 

After  data  was  taken  at  one  temperature  level,  the  power  to  the  main  heater  was 
increased  followed  by  equilibration  of  the  guard  cylinder  temperature  distribution.  When  a 
full  curve  of  contact  conductance  as  a  function  of  mean  interface  temperature  had  been 
generated  at  a  given  load  level  the  system  was  cooled.  In  some  cases,  data  points  were 
taken  during  the  cooling  cycle.  After  increasing  the  compressive  load,  another  contact 
resistance  curve  was  obtained.  Typical  variations  of  the  measured  quantities  are  given  in 
Table  V  along  with  the  resulting  uncertainty  in  the  contact  resistance  calculation. 

Several  important  points  can  be  distinguished  in  Table  V.  In  operation  of  the  apparatus, 
an  increase  in  mean  interface  temperature  is  attained  by  increasing  the  main  heater 
temperature,  the  "cold-end"  temperature  remaining  essentially  constant.  Thus,  as  the 
temperature  level  of  the  interface  increases,  two  effects  are  occurring  simultaneously: 

(1)  the  overall  temperature  drop  down  the  column  is  increasing,  and  (2)  the  specimen  thermal 
conductivity  is  decreasing  rapidly.  Since  the  overall  heat  flux  is  proportional  to  the  product 
of  these  quantities,  it  is  noted  in  this  column  of  Table  V  that  the  increasing  gradient  effect 
predominates  at  lower  temperatures  (increasing  flux)  while  decreasing  specimen 
conductivity  predominates  at  the  higher  temperatures  (decreasing  flux).  The  uncertainty  in 
the  heat  flux  involves  several  factors:  (1)  the  uncertainty  in  the  temperature  gradient 
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reflected  in  the  second  column  of  Table  V,  (2)  the  specimen  thermal  conductivity  uncertainty 
amounting  to  ±2%  (Section  II.  C.  2),  and  (3)  heat  exchange  with  the  thermal  insulation  due  to 
small  radial  temperature  gradients. 

At  lower  temperatures  the  uncertainty  in  the  thermal  contact  resistance  is  due  mainly  to 
errors  in  the  gradient  determination.  These  net  gradients  are  small  for  two  reasons 
(1)  the  heat  flux  is  low,  and  (2)  the  thermal  conductivity  is  high.  At  high  temperatures  the 
exchange  of  heat  between  the  specimen  column  and  the  insulation  becomes  important  because 
the  insulation  thermal  conductivity  is  increasing  (Appendix  III),  while  the  magnitude  of  the 
specimen  column  flux  is  decreasing  due  to  the  drop  in  conductivity.  The  combined  influence 
effects  (Table  V)  on  the  contact  resistance  measurements  is  minimized  in  the  intermediate 
temperature  range. 

A  complete  tabulation  of  the  experimental  result  is  given  in  Appendix  IV.  Graphical 
presentation  of  the  results  along  with  brief  qualitative  interpretations  are  given  in  the  next 
section. 

C.  GRAPHICAL  PRESENTATION  AND  QUALITATIVE  INTERPRETATION  OF  RESULTS 

The  graphical  presentation  of  measurement  results  giving  the  contact  conductance  as  a 
function  of  mean  interface  temperature  are  presented  in  the  next  10  Figures,  Figures  17-26. 
The  first  five  are  for  the  40  /x  in/40  fj.  in.  interface  at  compressive  load  levels  ranging  from 
about  5  psia  to  about  500  peia.  The  second  group  of  five  is  for  the  10  p  in.  /10  p  in.  interface 
covering  about  the  same  sequence  of  compressive  load  levels.  Formulation  of  an  interfacial 
heat  transfer  model  and  correlation  of  the  experimental  results  are  given  in  Section  W. 

The  following  brief  discussion  merely  points  out  the  general  features  and  trends  in  the  data. 

In  comparing  the  contact  conductance  ordinate  on  Figures  17-21  for  the  40  in.  /40  in. 
interface,  it  is  noted  that  the  contact  conductance  level  increases  rapidly  as  the  compressive 
load  increases.  The  hc  values  range  from  about  20  B1  U/hr-ft^-*  F  at  a  compressive  load  of 
approximately  5  psia  to  about  250  at  a  load  of  511  peia.  A  similar  load  dependence  of  h^  for 
the  smoother  10  fiin. /10  ^.in.  interface  is  noted  in  Figures  22-26.  Because  the  contacting 
surfaces  are  much  smoother  in  this  se'  and  sequence  of  tests,  the  level  of  interfacial 
conductance  is  higher  by  a  factor  of  between  3  and  4  at  any  given  load.  The  range  of  hg  is 
from  about  60  to  1000  for  compressive  loads  ranging  from  about  5  to  500  peia.  The  scatter 
in  the  experimental  data  tends  to  decrease  as  the  load  level  increases  in  accordance  with  the 
observations  of  the  previous  section. 
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Figure  17.  Contact  Conductance  as  a  Function  of  Interface  Temperature- 
40  /xin./40Ai.  in.  Interface,  5.02  peia  Compressive  Load 


Figure  18.  Contact  Conductance  a a  a  Function  of  Interface  Temperature 
40  ft.tn./40  n In.  Interface,  65.2  pelt  Compreestve  Load 
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Figure  19.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

40  fx  in./ 40  /Ain.  Interface,  149.3  psia  Compressive  Load 
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Figure  20.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

40  u.ln./40  ula.  Interface,  368.6  psla  Compressive  Load 
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Figure  21.  Contact  Conductance  as  a  Function  of  Interface  Temperature: 

40  ^iru/40  in.  Interface,  511.5  psia  Compressive  Load 
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Figure  23.  Contact  Conductance  as  a  Function  of  Interface  Temperature 
10  fj.  in./lO  fj.  in.  Interface,  64.7  psia  Compressive  Load 
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Figure  24.  Contact  Conductance  as  a  Function  of  Interface  Temperature 
10  ^iin./l0^.  in.  Interface,  148.7  psia  Compressive  Load 
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The  data  points  given  as  ,fX's"  in  Figures  18,  24,  and  25  are  cooling  curve  results.  These 
points  indicate  that  there  is  no  significant  difference  in  conductance  values  on  cooling  from 
the  high  temperature  exposure. 

Consider  now  the  general  shape  of  the  curves  in  these  ten  figures;  it  would  be  expected 
that,  if  the  heat  were  transferred  across  the  interface  only  at  the  actual  solid-solid  contact 
areas,  the  contact  conductance  would  exhibit  the  same  temperature  dependence  as  the 
beryllium  oxide  thermal  conductivity,  that  is,  a  decrease  in  contact  conductance  following 
a  1/T  dependency.  All  the  hc  versus  T  curves  exhibit  this  characteristic  at  low  temperatures. 
In  all  cases,  however,  the  hc  values  deviate  in  a  positive  sense  from  this  dependency  as  the 
temperature  level  increases.  As  noted  in  the  previous  section,  the  heat  flux  in  the  specimen 
column  increases  with  mean  interface  temperature  as  does  the  interfacial  temperature  drop. 
Since  the  radiative  heat  transfer  across  the  interfacial  gap  area  varies  as  Tm3  A  Tc  it 
would  be  expected  that  this  heat  transfer  mdchanism  becomes  important  at  higher 
temperatures.  Between  300  and  160(f  F,  which  represents  a  threefold  increase  in  absolute 
temperature,  the  interfacial  temperature  drop  across  the  40  /a  in.  /  40  ft  in.  contact 
increases  by  a  factor  of  between  4.  5  and  7,  depending  on  the  compressive  load.  The  magnitude 
of  the  radiative  contribution  thus  increases  by  a  factor  of  20  or  more.  The  situation  is 
similar  for  the  10  ft  in. /10  /a  in.  case,  the  interfacial  temperature  drop  increasing  bv  a 
factor  of  between  6  and  7  over  the  same  temperature  interval. 

The  net  result  is  a  predominance  of  the  radiative  contribution  at  high  temperatures  which 
determines  the  shape  of  the  hc  versus  T  curve  in  this  range.  In  most  cases  these  curves 
shift  to  a  positive  temperature  dependency;  a  quantitative  assessment  of  this  postulation  is 
given  in  the  next  section. 

Considering  the  shapes  of  the  curves  giving  the  40  ft  in. /40  ft  in.  results,  Figures  17-21, 
the  general  influence  of  the  compressive  load  on  the  relative  contributions  of  conductive 
and  radiative  transport  across  the  interface  region  can  be  inferred.  As  the  load  level 
increases,  it  is  clear  that  the  magnitude  of  solid  conductance  contribution  is  increasing 
rapidly;  the  magnitude  of  the  low  temperature  hc  values  verifies  this  observation.  The 
radiative  contribution,  on  the  other  hand,  is  not  directly  affected  by  the  compressive  load. 
Because  the  solid  conduction  contribution  is  comparatively  large  at  the  higher  loads  the 
interfacial  temperature  drops  are  lower  and  the  radiative  contributions  to  the  effective 
contact  conductance  are  consequently  lower.  Thus,  as  seen  in  Figures  20  and  21  which  give 
the  higher  load  data  for  the  40ftin/40  ft  in.  interface,  the  upturn  in  the  hc  curves  noted 


56 


AFML-TR-69-1 


at  lower  loads  has  been  eliminated.  In  effect  the  increased  solid  conduction  prevents  the 
generation  of  large  ATC  values  and  the  radiative  contributions  are  consequently  reduced. 


To  obtain  an  experimental  estimate  of  the  effective  gap  width  across  the  contact  region, 
additional  contact  conductance  tests  were  run  in  argon.  The  results  for  the  40  ^tin. /40  fx  in. 
and  the  10/iin./  10  /i.in.  interfaces  are  given  in  Figures  27  and  28,  respectively.  By 
comparing  the  vacuum  environment  thermal  contact  conductance  values  with  the  contact 
conductance  measurements  made  in  argon  at  the  same  compressive  load  and  temperature 
level,  the  contribution  of  gas  phase  conduction  to  the  total  heat  transfer  across  the  gap  can 
be  calculated.  Knowing  this  value,  the  effective  gap  width  can  be  determined  according  to 
Equation  25, 


8  « 


*  argon 

^c^argon  ^vacuum 


(25) 


These  data  were  obtained  at  a  temperature  low  enough  to  minimize  radiation  transport  across 
the  gap.  The  values  of  8  calculated  in  the  above  manner  were  used  as  a  check  of  the 
effective  gap  width  calculated  independently  from  the  surface  roughnesB-waviness  profile 
data.  These  results  are  discussed  in  the  next  section. 
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Figure  27.  Vacuum  and  Inert  Gas  Environment  Contact  Resistance  Results 
40  ^.in./40  ^itn.  Interface 
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Figure  28.  Vacuum  and  Inert  Gas  Environment  Contact  Resistance  Results: 
10  fiiu./lO  fj. in.  Interface 
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SECTION  IV 

DISCUSSION  OF  RESULTS 


A.  FORMULATION  OF  THE  CONTACT  MODEL 
1.  General  Contact  Characteristics 

The  temperature  discontinuity  observed  across  an  interface  transferring  heat  results 
from  the  imperfect  nature  of  contact,  as  shown  schematically  in  Figure  29  (Reference  1). 


Figure  29.  Flux  Field  Distribution  at  an  Interface 


60 


AFML-TR-69-1 


As  the  interface  is  approached  the  flux  lines  tend  to  converge  to  the  direct  solid-to-solid 
contact  points,  since  this  flow  path  generally  offers  considerably  less  resistance  to  heat 
transfer  than  the  surrounding  void  areas.  This  situation  would  be  altered,  of  course,  if  the 
void  areas  contained  a  gas  which  had  a  high  conductivity  or  if  the  temperature  level  were  so 
high  that  significant  radiative  transport  developed  across  these  gap  areas. 

The  convergence  of  the  flux  lines  to  the  solid-solid  contact  areas  produces  a 
"constriction  resistance"  and  an  interfacial  temperature  drop  to  compensate  for  it.  This  can 
be  noted  in  Figure  29  where  the  equipotential  (isothermal)  surfaces  (T)  and  (2)  ,  which 
are  orthogonal  to  the  flux  field,  are  seen  to  curve  away  from  the  contact  plane  at  the 
constriction  points. 

Early  investigations  in  this  field  either  neglected  surface  characterization  definition  or 
lumped  the  effects  of  surface  finished  on  interfacial  heat  transfer  into  an  effective  surface 
roughness,  usually  expressed  in  a  simple  RMS  or  CIA  value.  Beginning  in  1957-1958 
(e.g. ,  References  20  and  21),  as  the  results  of  surface  finish  analyses  were  applied  to 
studies  of  interfacial  heat  transfer,  it  became  apparent  that  such  a  simple  surface  charac¬ 
terization  procedure  was  entirely  inadequate  in  correlating  heat  transfer  across  real  surfaces 
in  contact.  Not  only  is  microscopic  roughness,  resulting  from  conventional  forming  and 
finishing  operations,  present  but  also  roughness  on  a  larger  wavelength  scale  such  as  surface 
waviness  or,  in  the  extremes,  surface  curvature  (Reference  1).  If  an  interface  is  made  up 
of  surfaces  which  are  both  rough  and  wavy,  relatively  large  contact  regions  governed  by  the 
waviness  will  be  formed.  Each  of  these  areas  will  consist  of  a  number  of  smaller  solid- 
to-solid  contact  points,  the  size  and  distribution  of  which  are  governed  by  the  surface 
roughness.  In  Figure  29,  the  central  contact  region  illustrates  this  characteristic,  since  it 
is  not  only  formed  as  a  result  of  the  contact  of  undulations  of  larger  scale  but  consists  also  of 
many  smaller  individual  contacts. 

As  a  result  of  extensive  surface  characterization  studies  described  briefly  in  Section  III.  A. , 
the  beryllium  oxide  specimens  used  in  this  investigation  were  found  to  be  both  rough  and  wavy. 
The  waviness  patterns  were  in  the  form  of  concentric  rings;  the  roughness  pattern  was  of 
random  lay,  as  expected  from  the  finishing  procedure,  and  consisted  of  individual  asperities. 

The  treatments  of  roughness  and  waviness  contributions  to  the  net  contact  resistance  are 
given  in  the  next  two  sections.  Following  this  is  the  presentation  of  the  combined  model  used 
in  correlating  the  experimental  contact  resistance  data  obtained  in  this  investigation. 
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2.  Surface  Roughness 

It  is  usually  assumed  that  the  surface  roughness  of  the  contacting  specimens  is  such 
that  a  group  of  point  contacts  with  approximately  circular  cross  section  are  formed  when  the 
specimens  are  brought  together.  With  this  approximation  it  is  then  necessary  to  develop  an 
expression  for  the  constriction  resistance  of  a  single  circular  cross-sectional  contact  as  a 
function  of  the  contact  area  radius.  Since  this  relationship  is  known  and  the  contact  points 
are  all  in  parallel,  the  total  resistance  of  the  group  of  contacts  is  found  simply  by  summing 
reciprocals  of  each  individual  resistance.  The  individual  contact  points  are  generally  assumed 
to  have  the  geometry  shown  on  the  left  in  Figure  30.  The  contact  area  radius  is  a,  the  average 
gap  width  is  8  ,  and  the  radius  of  the  element  delivering  heat  to  the  solid-solid  contact  is  b. 

An  expression  for  the  constriction  resistance  of  this  geometry  with  no  flow  across  the  gap 
area  was  worked  out  many  years  ago  by  Kottler  (Reference  22)  for  the  equivalent  electrical 
case.  Many  derivations,  producing  essentially  the  same  result,  both  exact  and  numerical, 
have  been  given  for  the  heat  transfer  case  as  discussed  in  Reference  1.  The  briefest  and 
most  direct  solution  was  given  by  Clark  and  Powell  (Reference  23),  who  considered  the  case  of 
steady-state  heat  transfer  as  being  a  gap  between  two  semi-infinite  regions  (1  and  2)  joined  by 
a  circular  contact  region  of  radius  a. 


The  Laplace  equation,  V  2t=0,  in  axisymmetric  cylindrical  coordinate  form  was 
written  for  each  region  with  the  following  coupling  boundary  conditions  along  the  contact  plane. 
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The  solution  by  "Hankel  transformation"  for  region  1  is  given  in  Equation  26,  where 
evaluation  of  the  integration  constants  has  been  completed  under  the  above  boundary  conditions. 
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A  sketch  of  this  temperature  field  and  the  orthogonal  heat  flux  paths  is  given  at  the  right 
in  Figure  30.  The  constant  temperature  surface  at  the  distance  from  the.  interface  where 
the  temperature  field  is  no  longer  distorted  due  to  the  constriction  is  represented  by 
If  kj  -  k0»  the  temperature  distribution  on  either  side  of  the  interface  is  symmetrical,  as 

mi 

depicted  in  Figure  30.  The  quantity  of  heat  flowing  from  one  region  to  the  other  is  given  by 
Equation  27,  using  the  relationship  expressed  by  Equation  26,  so  that 
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Figure  30.  Idealized  Contact  Element  and  Associated  Heat  Flux  -  Temperature  Field 
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The  constriction  resistance  is  then  defined  f»s, 

B  (T.  >i  =  0  .  __L_ 

Kc  -  Q  "  2k(  Q 


(28) 


Equation  28  is  the  fundamental  relationship  for  the  thermal  resistance  induced  in  a  conducting 
element  by  a  coaxial  region  of  reduced  cross-sectional  area  where  there  is  no  heat  transfer 
across  the  gap  area. 


The  diameter  of  the  "feeder"  element,  b,  does  not  appear  (Equation  28).  This  would  be 
expected  as  long  as  b  is  somewhat  greater  than  a  since  it  is  the  constriction  area  of  radius  a 
which  predominates  in  producing  the  temperature  field  distortion.  The  situation  whore  a 
approaches  b  has  been  studied  by  Roess  (Reference  24)  and  Clausing  and  Chao  (Reference  25). 
The  result  is  the  application  of  a  constriction  alleviation  factor  g(/3)  to  correct  Equation  28 
as  follows, 


9  <ff> 
C  2  k,  a 


(29) 


where 


q(/9)=  I  -  I  40925/?  ♦-  0.29591  /?3  -  + 

Equation  29  was  employed  in  calculating  the  surface  roughness  contributions  to  the  total 
interfacial  contact  resistance  using  experimentally  measured  thermal  conductivities  for  the 
beryllium  oxide,  kj,  and  values  of  a  and  (3  -  b/a  obtained  from  analysis  of  the  surface 
roughness  traces. 

3.  Surface  Wavincss 

The  configurations  of  the  contact  areas  between  the  beryllium  oxide  specimens 
produced  by  the  surface  waviness  were  in  the  form  of  concentric  rings.  The  contact  ring 
area  and  the  cross  section  shown  is  in  two-dimensional  Cartesian  coordinates.  The  contact 
parameters  a,  b,  and  8  were  determined  from  anal} sis  of  the  surface  waviness  profiles. 
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The  simple  analysis  of  the  previous  section  is  not  applicable  in  treating  surface 
waviness  contributions  to  the  total  contact  conductance  because  heat  transport  across  the  gap 
area  (a  <  r  <  b)  by  radiation  is  likely  to  be  important  at  higher  temperatures,  as  discussed 
qualitatively  in  Section  III.  C.  Radiation  contributions  were  neglected  in  the  case  of  the 
roughness  contact  resistance  calculation  for  two  reasons: 

(1)  The  temperature  field  distortions  caused  by  the  waviness  patterns  result  in  low 
temperature  drops  in  the  region  of  the  roughness  contacts,  while  the  temperature  drops 
across  the  gap  regions  between  waviness  contact  areas  are  relatively  high. 

(2)  The  percentage  of  the  total  interfacial  regions  occupied  by  roughness  contacts  is 
very  small  and  thus  the  radiating  surface  area  is  small.  The  total  radiation  contribution  is, 
of  course,  proportional  to  the  product  of  these  two  factors. 

To  enture  a  significant  degree  of  flexibility  in  handling  the  potentially  nonlinear  effects  of 
radiative  transport  in  the  gap  regions  between  waviness  contact  areas,  the  temperature  field 
in  the  contact  element  was  determined  numerically  by  the  finite  difference  solution. 

If  a  symmetrical  conducting  element  has  a  constriction  region,  the  centerline  of  the 
element  becomes  a  line  of  symmetry.  If  it  is  also  assumed  that  the  temperature  field  is 
symmetrical  with  respect  to  the  constriction  plane,  which  is  valid  when  kj  -  k 2,  then  study 
of  the  temperature  field  in  the  element  can  be  confined  to  a  quadrant  region  3uch  as  that 
sketched  in  Figure  31;  heat  flux  paths  are  indicated. 

After  the  steady-state  temperature  field  in  the  region  is  determined  with  the  boundary 
conditions  observed  experimentally,  the  effective  contact  conductance  can  be  inferred  from 
the  extent  of  distortion  of  the  heat  flux  path.  Radiation  boundary  conditions  occur  along  the 
surface 

c 

a  <  r  <  b  ,  0  <  z  <  -y 

A  typical  finite  element  network  is  shown  in  Figure  32,  In  this  Figure  there  are  72  nodal 
points,  whereas  the  networks  used  in  the  seryllium  oxide  analysis  averaged  about  500  nodal 
points.  All  the  different  types  of  nodal  po  nts  and  their  relative  locations  are  shown  in  this 
figure,  however.  There  are  six  types  of  nodal  points  including  fixed  temperatures  along  the 
surfaces  I,  J  ■  1  and  8. 
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Figure  31.  Contact  Element  Quadrant 

Two  approaches  can  be  used  in  formulating  the  nodal  equations  (the  results  are  identical). 
In  one  case,  the  finite  difference  form  of  the  equation 

V2  T  ;  0 

be  employed,  writing  equations  for  the  nodal  temperatures  in  terms  of  the  temperatures 
adjacent,  such  that  the  above  equation  is  obeyed.  In  the  other  case,  the  heat  fluxes  into  each 
node  can  be  summed  and  nodal  temperature  equations  derived  based  on  the  steady-state 
condition 


The  latter  approach  was  used  here  because  it  was  convenient  to  express  the  radiation  boundary 
conditions  in  terms  of  heat  flux  rather  than  temperature  equations. 

The  derivations  of  these  nodal  equations  are  given  in  Appendix  B.  The  results  are  given 
below  in  Table  VI.  The  equations  given  here  have  been  simplified  from  those  derived  in  the 
Appendix  by  assuming  that  Ax  =  At  -  t-  ,  throughout  the  nodal  network.  The  evaluation  of 
the  geometrical  and  emissivitv  shape  factors  in  the  heat  transfer  coefficient  equation  is  discussed 
in  the  next  section. 
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figure  32.  Finite  Element  Nodal  Point  Definitions 
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TABLE  VI 


NODAL  EQUATION  SUMMARY 


Type 

Nodal  Equation 

Internal 

T(I.J)=  ^[tII.J+D+TIKI.J-IM-KI  -+  l,J)+T  (I  —  1  ,  J  )  ] 

Line  of  Symmetry 
f  *  C.  Boundary 
\  -  •  Insulated 

T*I.  J>=  4“  [T  ( I  *  j  +  I)  +  T<I.J-I)+2T(I±  i.  j  )] 

Internal  Corner 

f  £ 
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It  is  noted  in  the  equations  of  Table  VI  that  the  parameter  describing  the  influence  of 
boundary  radiative  transport  on  the  temperature  field  is  the  dimensionless  group  h^^Vk. 

This  group,  of  Biot  modulus  form,  gives  the  ratio  of  boundary  radiative  transport  to  in-depth 
thermal  conduction  hR/(k/J).  For  a  given  overall  temperature  drop  in  the  element,  the 
magnitude  of  the  temperature  gradient  is  inversely  proportional  to  /  ,  hence  this  parameter 
is  grouped  in  the  conduction  term  (k/J). 

The  method  of  evaluating  the  finite  element  equations  consisted  of  assuming  an  initial 
temperature  distribution,  calculating  new  temperatures  for  each  nodal  point,  and  correcting 
the  original  temperature  by  a  small  percentage  of  the  difference  between  the  original  and 
newly  calculated  temperature.  This  is  basically  a  relaxation  method  with  the  use  of  a 
damping  factor  to  ensure  smooth  convergence  to  the  correct  temperature  field  as  repeated 
iterations  of  the  full  temperature  distribution  are  performed.  Convergence  stability  is  a 
concern,  of  course,  due  to  the  nonlinear  nature  of  the  radiative  boundary  condition;  the 
radiative  heat  transfer  coefficient  is  a  function  of  temperatures  which  vary  from  iteration  to 
iteration.  A  program  for  the  IBM  7094  computer  was  written  to  perform  the  iterative 
calculations.  Depending  on  the  number  of  nodal  points  used,  convergence  was  achieved  after 
100  to  600  iterations  of  the  entire  network.  Computer  output  consisted  of  the  final  temperature 
distribution  and  a  plot  of  the  nodal  residuals  throughout  the  network. 

As  an  illustration  of  the  procedure  used  in  calculating  the  contact  resistance  using  the 
finite  element  analysis,  the  results  of  two  cases  will  be  given  in  which  the  simple  network 
depicted  in  Figure  32  will  be  used.  In  the  first  case  the  temperature  along  the  surface  I, 

J  ■  1,  was  assumed  to  be  50(F  F,  while  that  along  the  surface  I,  J  ■  8,  was  assumed  to  be 
1500*  F.  Other  numerical  information  is  given  under  "Numerical  Evaluation"  in  Appendix  II. 
Nodal  residuals  for  the  entire  network  were  printed  after  2,  10,  50,  and  100  iterations. 
Maximum  values  of  these  residuals  were  8.  3,  3.  2,  0.  59,  and  0. 12%  of  the  mean  element 
temperature,  respectively.  The  temperature  field  after  100  iterations  is  shown  in  Figure  33. 
The  total  heat  flux  flowing  in  the  element  is  determined  using  Equation  30  and  the  temperature 
distribution  away  from  the  interface  (e.g. ,  the  1400  and  150(f  F  isotherms)  is 

"total  t  AI  (-%-)T  (30) 

where  A  T  is  10(f  F  in  this  case  and  A  Y  is  the  distance  between  these  two  isotherms. 

The  radiative  heat  transfer  across  the  gap  area  was  calculated  similarly  using  Equation  31 
and  the  temperature  distribution  along  the  gap  surface. 
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radiation  .  . 

1=4 


(T i.j>  Ai  ATi 


(31) 


The  heat  flow  across  the  solid  contact  area  is  calculated  as  the  difference  between  the 
above  quantities  or  can  be  determined  directly  by  using  the  temperature  field  data  in  the 
vicinity  of  the  solid  contact  area  along  with  an  equation  of  the  same  form  as  Equation  30. 

For  this  case  only  3%  of  the  heat  flows  across  the  radiation  gap  which  produces  the  large 
distortion  of  the  temperature  field  noted  in  Figure  33. 


To  test  the  stability  of  the  relaxation  iteration  technique  the  influence  of  radiative  heat 
transport  in  the  gap  was  accentuated  by  arbitrarily  increasing  the  constants  in  the  equation  for 
(  ^Radiation)  ^y  a  ^actor  °f  1000*  This  artifically  high  effect  of  radiative  transport 
produces  the  temperature  field  distribution  shown  in  Figure  34.  After  100  iterations  the 
maximum  nodal  residuals  in  this  case  were  on  the  order  of  0.0001%.  Thus,  excellent 
stability  is  attained  with  the  damping  factor  ever  in  the  presence  of  large  radiation  contributions. 
The  form  of  the  temperature  field  in  Figure  34  indicates  that  heat  is  preferentially 
transferred  across  the  gap  region.  Indeed,  heat  flux  calculations  using  Equations  30  and  31 
indicated  that  65%  of  the  heat  flowing  in  the  element  was  transferred  across  this  region. 

The  Biot  modulus,  ^/k.  for  these  two  cases  differ  by  a  factor  of  greater  than  400. 

For  the  temperature  distribution  given  in  Figure  33  the  value  was  0.  Oil;  for  the  distribution 
in  Figure  34  the  value  was  4.66. 


To  calculate  the  thermal  contact  conductance  from  these  temperature  field  distributions 
the  following  procedure  was  used.  First,  the  temperature  drop  required  to  transfer  the  same 
amount  of  heat  in  an  element  with  no  discontinuity  was  calculated.  For  the  distribution  shown 
in  Figure  33  the  total  heat  flowing  (with  the  constriction  present)  was  4,243  BTU/hr.  The 
overall  A  T  as  shown  in  the  Figure  was  1500  -  500  ■  1000*  F.  The  temperature  drop 
required  to  transfer  this  amount  of  heat  in  this  element  with  no  constriction  was  668*  F. 

Thus,  the  temperature  drop  developed  due  to  the  presence  of  the  contact  was  ATC  =  1000  -  668 
•  3321*  F.  Stated  another  way,  the  presence  of  the  constriction  causes  displacement  of  the 
1168*  F  isotherm  (1500  -  332  »  1168*  F),  shown  as  the  dashed  curve  in  Figure  33,  to  the 
I,  J  =  8  plane.  TTiis  A  TJ,  divided  by  the  heat  flux  crossing  the  solid  and  gap  areas  gives 
the  thermal  contact  resistance  of  each  region.  The  total  contact  conductance  for  the 
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Figure  33.  Isotherm  Distribution  in  the  Contact  Element  Quadrant  with  Normal 
Radiation  Heat  Transfer  Coefficient 
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distribution  given  in  Figure  33  was  (^c)^  -  150. 93  BTU/hr-ft^-*  F.  The  value  for  the 
distribution  in  Figure  34  was  (h)j  =  -466.7.  The  negative  quantity  is  expected  since  the 
artifically  enhanced  gap  transport  is  a  higher  conduction  region  than  the  element  material 
itself  (hj^i’/^  >  1). 

The  application  of  a  finite  element  analysis  to  the  constriction  resistance  problem  has 
not  been  presented  in  the  literature.  Cetinkale  and  Fishenden  (Reference  26)  indicate  that 
the  Southwell  relaxation  method  was  used  to  obtain  contact  resistance  values  in  a  two- 
dimensional  cylindrical  contact  element.  Their  analysis  probably  involved  the  use  of  an 
extensive  grid  network  although  the  details  are  not  discussed.  These  authors  elaborated  in 
much  greater  detail  on  the  use  of  an  exact  analytical  method  of  approximating  the  temperature 
field  in  the  proximity  of  a  constriction  and  the  explicit  formulation  of  contact  resistance 
equations  based  on  the  model. 

4.  The  Combined  Model 

Analysis  of  the  surface  roughness  and  waviness  profiles  indicated  that  the  beryllium 
oxide  test  surfaces  exhibited  significant  degrees  of  both  of  these  forms  of  surface  irregularity. 
The  interfacial  heat  transfer  model  included  both  of  these  effects.  The  first  experimental 
evidence  that  more  than  the  surface  roughness  contributed  to  the  mean  interfacial  gap 
(of  thickness  8  )  came  from  contact  resistance  measurements  in  argon.  Equation  25  and  the 
experimental  results  presented  in  Figures  27  and  28  for  the  40  and  10  ^.in.  interface 
respectively  were  used  to  calculate  an  effective  gap  width.  The  average  values  were  119.6 
and  756. 3  in.  for  the  10  and  40  ^tin.  CLA  interfaces,  respectively.  Where  only  surface 
roughness  contributions  are  concerned,  the  effective  gap  widths  calculated  from  the  surface 
profiles  were  about  50  ft  in.  for  the  smoother  interface  and  100  jj.  in.  for  the  rougher  40  yu.in. 
interface.  The  disparity  between  these  results  suggested  the  presence  of  surface  waviness 
patterns  which  would  tend  to  increase  considerably  the  effective  8  .  Estimated  waviness 
contributions  are  shown  in  Table  VII  along  with  the  other  results  mentioned  above. 

The  waviness  contributions  add  another  60-70%  to  the  effective  8  although  the  values 
still  do  not  equal  the  estimates  obtained  from  the  contact  resistance  data  in  argon.  The  mean 
free  path  of  the  argon  at  the  temperature  and  gas  phase  pressure  of  these  tests  was  5. 39  fi  in. 
the  mechanism  of  gas  phase  conduction  should  then,  be  free  of  Knudsen  transport  regime 
effects.  It  would  be  expected  that  the  accommodation  coefficient  between  the  relatively  inert 
beryllium  oxide  and  argon  would  be  less  than  unity.  The  apparent  increase  in  the  effective 
gap  width  due  to  this  characteristic  has  been  discussed  by  Kennard  (Reference  27)  and  by 
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TABLE  VII 

EFFECTIVE  GAP  WIDTH  COMPARISONS 


8  ,  /tin. 

Interface 

Roughness 
Contribution 
per  Specimen, 

&R 

Waviness 

Contribution 
per  Specimen, 

8  w 

Total 

<28r±2Sw) 

Measured  in 

Argon  using 
Equation  25 

10  /tin. 

25 

15 

80 

119.6 

40  ft  in. 

120 

80 

400 

756.3 

Rapier,  Jones,  and  McIntosh  (Reference  28).  The  apparent  increase  in  the  effective  gap 
width  due  to  accommodation  coefficients  less  than  unity  is  termed  the  accommodation  jump 
distance,  g,  by  these  authors.  It  is  proportional  to  the  molecular  mean  free  path  and  the 
quantity  j(2-  /t  )//t  j  where  /t  is  the  accommodation  coefficient.  The  value  of  g, 
in  microinches,  which  would  add  directly  to  the  "Total"  column  ir.  Table  VII,  ranged  between 
28.  0  and  61.  0  depending  on  the  value  of  ft  ;  values  from  0.  5  to  0. 85  were  considered. 


These  results  along  with  the  surface  profile  traces  clearly  indicate  the  importance  of 

both  roughness  and  waviness  effects  for  the  beryllium  oxide  surfaces  studied  here.  The 

quantitative  combined  model  used  in  generating  ar  '.lytically  predicted  curves  of  contact 

conductance  as  a  function  of  temperature  (and  r^mpreasive  load)  is  shown  in  Figure  35.  The 

contact  rings  i  and  J  produced  by  the  waviness  pat*  ’•ns  each  consist  of  a  number  of  contact 

points  produced  by  the  roughness  patterns.  The  the  -nal  resistances  induced  by  the  waviness, 

(Rw,  Rw  )j  and  (Rw,  R\y)i»  are  in  serles  with  those  induced  by  the  roughness,  (Rr{rh 
I  II  I  II  f  ’ 


^Rr^hh  and  ^Rr^h  *  Rr^h^  j  »  whlle  the  radiation  gap  resistances  R^  and  R^d  are  in 
parallel  with  both.  The  total  thermal  contact  resistance  of  the  network  shown  in  Figure  35  is 


given  by  Equation  32. 
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<VT 


U) 


[RRad+RRod]+[lRw+R"+Rrgh  4Rrgh,j  +<Rw4Rw  4 Rrgh+Rrgh*.  ] 


(32) 


Ui 


where 

£  =  RRod  RRod 

»  i 

"  =  <Rw4Rw  4Rfgh4Rrgh)i  (Rw  +  Rw4Rrgh+Rroh»( 

i  d  i  n.  i  n  i  u 

The  waviness  and  radiative  thermal  contact  resistances  R^  and  were  determined 

by  the  methods  outlines  in  Section  III.  The  roughness  thermal  contact  resistances,  Rrgh> 

were  calculated  using  the  methods  outlined  in  Section  II.  In  Section  V  comparisons  are  made 
of  these  analytically  predicted  contact  resistances  with  the  experimentally  measured  contact 
resistance  versus  mean  interface  temperature  results. 

B.  COMPARISON  OF  ANALYTICALLY  PREDICTED  RESULTS  WITH  EXPERIMENTAL 
MEASUREMENTS 

To  generate  an  analytically  predicted  (hc)t  versus  mean  interface  temperature  (Tm) 
curve,  three  types  of  input  data  are  required 

a.  Surface  geometry  data  (roughness  and  waviness) 

b.  Temperature  dependent  physical  property  data  (thermal  conductivity,  emittance, 
and  compressive  strength)  for  the  beryllium  oxide. 

c.  Gross  temperature  distribution  data  Tm  and  an  estimate  of  the  range  of 
temperatures  found  in  the  vicinity  of  the  contact  for  a  given  Tm. 

1.  Input  Parameters 

Analysis  of  the  surface  profile  data  produced  the  set  of  geometrical  parameters  for 
the  10 and  40  in  surfaces  summarized  in  Table  VII. 

The  thermal  conductivity  data  used  was  measured  on  the  test  specimens  and 
reported  in  Figure  10.  The  compressive  strength  data  for  BeO  in  the  density  and  purity 
range  used  here  was  taken  from  the  recent  ceramic  materials  compendium  issued  by 
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Battelle  Memorial  Institute  (Reference  29).  Emittance  data  on  beryllium  oxide  was  run  during 
the  course  of  this  investigation  on  specimens  in  the  same  surface  roughness  range  as  the  contact 
specimens.  Spectral  normal  emittance  was  determined  using  a  Perkin-Elmer  spectrophotometer 
in  accordance  with  the  methods  described  in  Reference  30.  The  numerical  integration  required 
in  the  conversion  of  this  spectral  data  to  the  total  hemispherical  emittance  values  needed  in 
the  heat  transfer  calculation  are  discussed  in  Appendix  m. 

The  radiative  heat  transfer  coefficient  given  in  Table  VI  requires  evaluation  of  the 
geometrical  and  emittance  shape  factors,  and  ,  respectively.  From  Table  V1IT  it 
is  seen  that  the  gap  thickness  to  gap  width  ratio  varies  from  1.  OxlO-4  to  2.  7xl0_3;  thus,  the 

assumption  that  radiative  transport  occurs  between  infinite  parallel  planes  is  a  good  approximation 
(Reference  31).  This  leads  to  the  conclusion  that  J^0  =1.0  and  that, 


<  t 

I  2 


as  long  as  the  radiative  transport  is  in  the  geometrical  optics  regime.  The  longest  photon 
wavelength  of  interest  in  ti  e  Planck  radiation  distribution  occurs  at  low  interface 
temperatures  and  was  on  the  order  of  10  fx  in.  Thus,  the  assumption  of  geometrical  optics, 
on  which  the  derivation  of  Equation  33  is  based,  is  valid. 

Evaluation  of  Equation  33  is  performed  internally  during  the  execution  of  the  IBM  7094 
finite  element  computer  program  using  the  varying  radiation  node  temperatures  and  tabular 
«TH*f(T)  data  in  core  storage  to  determine  appropriate  «,  and  «2  values. 

2.  Analytically  Developed  Contact  Conductances 

The  waviness  pattern  dimensions  in  Table  VIII  establish  the  basic  geometry  of  the 
contact  element  to  be  analyzed  in  the  finite  element  computer  program.  The  spacing  of  the 
nodal  points  as  determined  by  the  parameter  J!  was  chosen  small  enough  to  give 
numerous  temperature  points  in  regions  of  maximum  temperature  field  distortion.  The  size 
of  the  nodal  networks  varied  from  20x28  for  the  10  fj.  in.  interface  to  20x38  for  the  40  fJ.  in. 
interface.  Values  of  2  varied  from  0.0039  to  0.0044  inches. 

With  the  geometry  of  the  network  established,  it  was  necessary  to  aseign  boundary 
temperatures  along  the  fived -temperature  nodes.  The  relative  position  of  these  nodes  is 
given  in  Figure  32.  The  assignment  of  these  temperatures  is  not  particularly  critical  at  low 
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temperatures  as  long  as  the  mean  contact  element  temperature  near  the  interface  is  approximately 
equal  to  T  ,  the  experimentally  observed  mean  interface  temperature.  This  is  because  the 
contact  conductance  is  nearly  heat  flux  independent  at  low  temperatures;  variation  of  the 
fixed-node  temperatures  at  constant  Tm  changes  only  the  net  flux  in  the  contact  element.  At 
high  temperatures  the  contact  conductance  becomes  temperature  level  and  temperature 
gradient  dependent  due  to  the  radiative  contributions;  thus,  care  is  required  in  assigning 
boundary  temperatures  to  produce  Tm  and  temperature  gradient  values  which  correspond  to 
those  achieved  experimentally. 

To  check  the  computer  execution  of  the  finite  element  program  a  hand  calculation  of  one 
iteration  for  the  entire  network  was  performed.  The  nodal  temperatures  and  the  nodal 
residuals  agreed  exactly  with  those  generated  in  the  computer  calculation.  During  the  course 
of  the  successive  iterations,  selected  nodal  indexes,  iteration  indexes,  nodal  residuals, 
nodal  temperatures,  and  heat  transfer  coefficients  were  printed  out.  This  served  as  a 
check  of  calculation  progression  and  allowed  manual  calculation  of  the  heat  transfer 
coefficients  for  comparison  with  computer-calculated  values.  In  all  cases  agreement  was 
exact.  For  the  larger  20X30  node  fields,  convergence  was  fairly  slow  and,  thus,  printout 
of  the  complete  nodal  residual  network  was  made  to  determine  the  number  of  iterations 
required.  On  the  average,  600  iterations  were  sufficient  to  reduce  the  maximum  residual 
to  no  more  -han  0.  002%  of  the  average  temperature  of  the  element. 

To  illustrate  the  application  of  the  combined  contact  model,  a  set  of  (hc)-p  versus  T 
curves  for  the  10  and  40  in.  interfaces  at  a  given  load  level  were  analytically 
generated  for  comparison  with  experimental  measurements  at  the  same  load  level.  The 
low-load  level  curves  for  each  interface,  given  in  Figures  17  and  22  were  chosen.  In  the 
case  of  the  analytical  prediction  for  the  40  ft  in.  interface,  it  was  assumed  that  the 
contact  area  increased  through  brittle  fracture  at  the  peaks  until  the  load  of  3.  94  lbs  was 
supported.  For  the  10  fx  in.  interface  it  was  assumed  that  the  average  load  supported  by  the 
asperities  was  less  by  a  factor  of  6.  2,  the  ratio  of  the  contact  areas  versus  displacement 
curves,  for  the  two  surfaces  (Figure  16). 

This  assumption  is  based  on  the  conclusion  that  very  few  of  the  10  u  in.  interface  peaks 
undergo  brittle  fracture  in  supporting  the  fairly  low  3.68  lbs  load  applied  in  this  test.  The 
surface  roughness  profiles  (Figure  14)  and  the  asperity  geometry  (Table  VIII)  for  this 
interface  indicate  that  the  polishing  employed  produced  contact  areas  with  comparatively 
broad  peaks.  The  results  of  the  roughness  and  waviness  contact  resistance  calculations  are 
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summarized  in  Table  tX.  These  values  were  employed  in  calculating  the  total  interfacial 
contact  resistance  and  its  reciprocal,  the  total  interfacial  contact  conductance  (hc)-p,  using 
Equation  32.  The  values  were  based  on  total  heat  transfered,  Q  =  BTU/hr,  rather  than  on 
heat  flux,  q  *  BTU/hr-ft  .  After  completion  of  the  calculations  the  values  were  converted 
to  a  flux  per  unit  area  basis.  In  Table  EX  there  is  a  consistent  trend  of  increasing  roughness 
and  waviness  contact  resistance  values  as  the  temperature  increases.  This  is  due  mainly  to 
the  decreasing  thermal  conductivity  of  the  beryllium  oxide.  The  roughness  and  waviness 
contact  resistances  are  about  the  same  order  of  magnitude.  Also  consistently  observed  is  a 
decrease  in  the  radiative  contact  resistance  with  increasing  temperature.  The  decrease  is  by 
a  factor  of  between  15  and  20  over  the  temperature  range  studied.  Such  a  decrease  in 

_o 

resistance  is  expected  due  to  the  T  dependence  of  this  term. 

The  final  results  of  the  analytical  calculation  using  Equation  32  are  compared  with  the 
experimental  measurements  on  the  10  and  40  fi  in.  interfaces  in  Figures  36  and  37, 
respectively.  The  predicted  conductive  and  radiative  contributions  to  the  total  interfacial 
conductance  are  also  indicated.  As  inferred  qualitatively  in  Section  III.  C. ,  radiative 
contributions  are  much  more  important  for  the  rougher  40  fx  in.  interface,  the  conductive 
and  radiative  contributions  becoming  equal  at  about  95(f  F  for  this  interface;  for  the  smoother 
interface,  the  radiative  contribution  does  not  equal  the  conductive  contribution  until  a 
temperature  of  145 (T  F  is  reached.  In  absolute  magnitude,  the  radiative  contributions  are 
about  the  same  for  the  two  interfaces.  However,  the  solid  contribution  is  4.  5  times  greater 
for  the  smooth  interface,  thus  decreasing  the  relative  importance  of  radiative  transport. 

'rhe  dimensionless  grouping  of  Biot  modulus  form,  hR  j? / k,  introduced  earlier,  was 
calculated  for  each  of  the  analytically  predicted  data  points.  Calculation  of  this  group  termed 
the  Interfacial  Heat  Transport  Modulus,  Hj,  is  summarized  in  Table  X.  The  characteristic 
length,  1,  is  the  horizontal  distance  between  waviness  contact  areas,  which  is  a  quantitative 
measure  of  the  distance  from  the  interface,  in  which  there  is  significant  perturbation  of  the 
temperature  field. 

3 

In  general,  h^  -  f(Tm  )  and  k  -  f(l/Tm),  thus  it  would  be  expected  that  the  interfacial 
heat  transport  modulus  should  exhibit  a  T*  functionality, 

HI  =  hR^  /  k  =  f  (T^  )  (3  4) 
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Figure  36.  Experimental-Analytical  Contact  Conductance  Comparisona:  10  /i  in.  Interface 
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Figure  37.  Experimental- Analytical  Contact  Conductance  Comparisons:  40/xtn.  In 
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where  Tm  is  in  degrees  Rankine.  The  values  in  Table  X  are  plotted  on  log-log  coordinates  in 
Figure  38.  The  Hj  values  clearly  define  straight  lines.  The  slopes  are  the  same  and  are  just 
slightly  above  4.  0,  in  both  cases  equaling  4.4.  That  the  curves  for  the  two  interfaces  do  not 
coincide  but  are  offset  by  a  constant  amount  is  expected.  This  is  because  the  absolute 
magnitude  of  the  solid  conduction  contributions  at  low  temperature  are  different  for  each 
interface: 

k.oMin.=  ‘V’m'  <V  U/T„,)  ,  C,  *  C, 

The  Hj  value  is  useful  because  it  gives  a  direct  measure  of  the  relative  contributions  of 
conductive  and  radiative  heat  transport  across  the  interface.  Knowing  the  value  for  a  given 
interface  at  a  given  temperature  allows  estimation  of  the  value  at  other  temperatures  through 
the  (Tm)4.4  dependence. 

Clausing  and  Chao  (Reference  25)  have  expressed  the  interfacial  contact  conductance  in 
dimensionless  form  using  the  ratio  of  (hc)T  times  a  characteristic  surface  dimension  to  the 
mean  thermal  conductivity  of  the  contacting  specimens.  These  authors  term  this  ratio  a  Biot 
number  although  it  has  an  entirely  different  physical  significance  from  the  interfacial  heat 
transport  modulus,  Hj.  Clausing  and  Chao  considered  the  influence  of  roughness  and 
waviness  on  solid  conduction  across  the  interface.  Thus  the  Biot  number  they  employ  is 
basically  a  ratio  of  solid  conduction  across  the  interface  to  conduction  in  the  bulk.  The  Hj 
grouping  defined  here  is  used  in  the  more  conventional  sense  of  the  Biot  modulus  (Reference  32). 
That  is,  it  is  defined  as  the  ratio  of  a  surface  heat  transfer  coefficient  to  a  bulk-media 
transport  property. 

The  predicted  contact  conductance  curves  in  Figures  37  and  38  fall  an  average  of  11% 
below  the  measured  values  although  they  are  within  the  scatter  of  the  experimental  results 
which  was  estimated  to  be  ±12  to  ±23%  (Table  V).  The  slightly  lower  predicted  values  are 
attributed  to  th°  fact  that  the  contact  model  somewhat  underestimates  the  solid  conduction 
contribution  du  •  to  the  reduction  in  compressive  strength  of  the  beryllium  oxide  with 
increasing  temperature.  Because  of  the  high  radiative  contributions  at  elevated  temperatures, 
both  predicted  curves  are  sensitive  to  the  emittance  values  in  this  range.  The  possible 
uncertainties  in  this  data  at  elevated  temperature  is  indicated  in  Figure  43  of  Appendix  III. 

The  emittance  data  measured  during  this  investigation  falls  between  other  data  reported  in 
the  literature,  which  exhibits  a  very  considerable  spread. 
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The  analysis  used  here  and  the  results  obtained  suggested  a  simpler,  partially  empirical 
method  of  obtaining  predicted  curves  of  (hc)-p  as  a  function  of  Tm.  This  simplified  analysis 
is  given  in  the  next  section. 

C.  A  SIMPLIFIED  ANALYTICAL  PREDICTION  PROCEDURE 

During  the  experimental  measurements  a  value  of  ATc  is  obtained  directly  by  gradient 
extrapolation  to  the  interface.  There  certainly  are  temperature  variations  radially  along  the 
gap  surfaces  but  the  finite  element  analysis  showed  these  to  be  fairly  small.  Thus,  the 
measured  A  Tc  gives  a  representative  average.  The  radiative  transport  across  the 
interface  could  then  be  approximated  fairly  accurately  by  Equation  35, 

"Rod*  'S'  (4T-"’  ’AT‘  138 

and 

lhC  >Rod  “  I 

A  great  amount  of  surface  profile  analysis  is  required  to  arrive  at  an  independent  estimate 
of  the  interfacial  contact  area,  and  ultimately,  the  solid  conduction  contribution.  Now,  if 
the  radiative  contributions  to  the  total  contact  conductance  are  small  at  low  temperatures, 
it  can  usually  be  assumed  that  the  total  inter  facial  heat  transfer  is  due  only  to  conduction 
across  the  solid  contacts.  Knowing  the  beryllium  oxide  thermal  conductivity  temperature 
dependency  makes  it  reasonable  to  assume  that  the  change  in  the  solid  conduction  contribution 
to  the  total  contact  resistance  with  temperature  follows  the  same  dependency.  If  the  total 
low  temperature  contact  conductance  is  defined  as  (hc)0  then  the  solid  conduction 
contribution  as  a  function  of  temperature  can  be  approximated  as 

">«  ■  <»«>«  +  ">«!,,„  -  lhc„  1/Tm 

Using  the  experimentally  determined  values  of  T  ,  AT  and  (h  )  and  Equations  35 

m  c  c  o 

and  36,  the  predicted  contact  conductance  curves  shown  in  Figures  39  and  40  were  generated 
for  the  10  and  40  fi  in.  interfaces,  respectively. 

The  shapes  of  the  predicted  curves  in  these  figures  are  slightly  different  from  those 
obtained  from  the  exact  analysis  described  earlier,  but  still  the  general  similarity  and  the 
agreement  between  experimental  and  predicted  results  is  striking  considering  the  simplicity 
of  this  analysis.  As  stated  initially,  however,  an  experimental  contact  resistance  measurement, 
(hc)0,  is  required  to  use  the  analysis. 
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Figure  39.  Experimental-Analytical  Contact  Conductance  Comparisons. 
Simplified  Analysis  -  10^.  in.  Interface 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  section  conclusions  and  recommendations  are  based  on  the  extensive  measurements 
and  analyses  conducted  in  assessing  the  performance  of  the  experimental  equipment  and  in 
studying  the  character  of  heat  transport  across  beryllium  oxide  interfaces,  as  already 
discussed  in  this  report.  The  conclusions  and  recommendations  are  presented  bv  section  topic 
as  given  earlier. 

A.  CONCLUSIONS 

1.  Design,  Construction,  and  Qualification  Testing  of  Experimental  Equipment 

The  design  philosopy  of  providing  a  large  lateral  thickness  of  thermal  insulation 
for  the  relatively  high  conductivity  specimen  column  to  minimize  radial  losses  while  at  the 
same  time  maintaining  only  approximate  thermal  balancing  requirements  on  the  guard  heater 
was  proved  sound.  Continual  operation  with  hot  side  temperatures  of  over  2000*  F  was 
achieved.  Heat  losses  were  held  to  within  5%  with  guard  mismatching  ranging  from  -*-23  to 
-11%,  while  thermal  conductivity  data  was  obtained  in  excellent  agreement  with  data  reported 
in  the  literature. 

The  beryllium  oxide  thermal  conductivity  data  measured  with  the  equipment 
designed  and  built  during  this  program  had  a  precision  to  within  better  than  ±2%.  The 
temperature  dependence  of  the  results  was  in  exact  accord  with  theoretical  predictions  and 
the  literature  values  (bracketed  data)  available  on  specimens  in  the  same  purity  and  density 
range.  Estimated  accuracy  of  the  data  was  ±3%.  The  thermal  conductivity  data  run  on 
tantalum  as  an  additional  check  of  equipment  performance  exhibited  a  precision  of  ±1.5^. 
Available  literature  values  for  this  metal  exhibit  a  wide  scatter;  however,  the  data  obtained 
in  this  investigation  agreed  with  recommended  values  within  ±2%  up  to  750*  F,  and  deviated 
no  more  than  7.5%  at  120 (f  F. 

Analytically  predicted  heat  losses  in  the  apparatus  agreed  within  a  factor  of  1.  2  to  3 
with  experimental  measurements.  Uncertainties  in  rvailable  physical  property  data  and  in 
true  temperatures  at  the  insulation  boundaries  prevented  exact  application  of  these  analytical 
equations  which  could  easily  be  parameterized  to  produce  general  design  equations  for  the 
cut-bar  type  of  conductivity-contact  resistance  apparatus. 
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2.  Experimental  Generation  and  Analyses  of  BeO  Interfacial  Contact  Resistance 

Detailed  surface  characterization  is  required  which  includes  precise  analysis  of 
relatively  large  scale  surface  undulations  (waviness  and  curvature)  as  well  as  the  more 
conventionally  studied  microscopic  surface  irregularities  (roughness).  For  the  carefully 
prepared  interfaces  studied  here,  roughness  and  waviness  contribute  about  equally  in 
determining  the  nature  of  solid  phase  heat  transport  across  the  interfaces. 

The  physical  character  of  the  10  /i  in.  CLA  and  the  40  ^.in.  CLA  surfaces  studied 
here  were  fundamentally  different.  The  10  in.  specimens  had  relatively  flat  asperity 
surfaces  reflecting  the  polishing  procedures  used,  while  the  40  fx  in.  specimens  exhibited 
a  rougher  peak  pattern  due  to  the  vapor-blasting  operation  employed  in  surface  preparation. 
These  differences  in  geometrical  character  lead  to  significantly  different  interfacial  heat 
transfer  patterns,  which  were  consistently  explained. 

The  thermal  contact  conductance  results  exhibited  several  consistent  trends: 

a.  The  values  for  the  smoother  10  /u.in.  interface  were  higher  by  a  factor  of 
3  to  4  than  those  for  the  40  in.  interface  at  comparable  compressive  load  levels. 

b.  The  values  increased  rapidly  with  compressive  load. 

c.  Radiative  heat  transport  across  the  interfaces  was  important  at  high  temperatures, 
particularly  at  low  compressive  load  levels.  This  contribution  increases  by  a  factor  of  up 

to  20  in  the  temperature  range  from  200  to  1600*  F  and  was  a  major  contributor  to  the  net 
interfacial  heat  transport. 

d.  The  contact  conductance  curves  were  repeatable  on  cooling.  And, 

e.  The  compressive  load  on  the  interface  serves  to  reduce  the  relative  contributions 
of  radiative  transport  indirectly  by  preventing  the  generation  of  large  A  Tc  values  at 
higher  temperatures. 

An  estimate  of  the  effective  interfacial  gap  width  was  obtained  through  thermal  contact 
resistance  measurements  in  argon.  These  data  confirmed,  by  direct  measurement,  the 
importance  of  surface  waviness  on  the  pattern  of  interfacial  heat  transport. 

A  combined  roughness -waviness  interfacial  heat  transfer  model  was  formulated  which 
successfully  correlated  the  experimental  measurements  within  10%,  including  the  contributions 
of  high  temperature  radiative  transport  across  the  contact.  The  thermal  resistances  induced 


91 


AFML-TR-69-1 


by  surface  waviness  were  in  series  with  those  induced  by  the  roughness,  the  combination 
being  in  parallel  with  the  thermal  resistance  across  the  radiation  gap  regions. 

The  nonlinear  character  of  the  radiative  transport  effects  were  successfully  incorporated 
in  a  finite -element  analysis  that  exhibited  stable  convergence  over  wide  ranges  of  radiation 
heat  transfer  coefficient. 

The  relative  contributions  of  radiative  and  conductive  transport  across  the  interface  and 
their  combined  influence  on  the  temperature  field  distortion  were  characterized  bv  a 
dimensionless  interfacial  heat  transport  modulus,  Hj  *  hj^^/k.  This  group  was  proportional 
to  the  4.4  power  of  the  mean  interface  temperature,  Tm,  expressed  in  °  R;  a  4.  0  dependence 
was  predicted  theoretically  for  the  beryllium  oxide  interface. 

Analyses  using  the  combined  contact  heat  transport  model  showed  quantitatively  that 
radiative  contributions  were  much  more  important  for  the  rougher  40  /x  in.  interface, 
equalling  the  conductive  contribution  at  950°  F  for  the  low  compressive  load  case.  For  tlx* 

10  ft  in.  interface,  these  contributions  were  not  equal  until  an  interface  temperature  of 
145(f  F  was  reached. 

A  simplified  analytical  procedure  for  predicting  the  temperature  dependency  of  the  thermal 
contact  conductance  was  applied  successfully,  giving  essentially  the  same  results  at  the  more 
detailed  combined  model  analysis.  This  simplified  procedure  was  empirical  in  nature, 
however,  requiring  one  contact  resistance  measurement  for  a  given  interface  at  a  given 
load  level. 

B.  RECOMMENDATIONS 

1.  Experimental  Equipment 

The  apparatus  should  be  modified  to  include  a  secondarv  heater  between  the  bottom 
of  the  specimen  column  and  the  water-cooled  copper  pedestal  (Figure  2).  This  modification 
would  allow  independent  variation  of  the  mean  interface  temperature  and  the  heat  flux  in  the 
specimens.  Operation  in  this  fashion  is  not  possible  unless  there  is  independent  control  of  the 
temperature  at  the  bottom  of  the  specimen  column. 
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The  analytical  heat  loss  equations  developed  in  Section  II.  B.  3  should  be  used  to  generate 
generalized  apparatus  design  curves  such  as  the  following: 

a.  qlogg  as  a  function  of  insulation  thermal  conductivity  with  specimen  conductivity 
as  a  parameter. 

b.  q]oss  as  a  function  of  insulation  to  specimen  diameter  ratio,  with  specimen  to 
insulation  conductivity  ratio  as  a  parameter.  And, 

c.  q^  as  a  function  of  guard-specimen  column  temperature  gradient  mismatch 
with  both  of  the  parameterizations  mentioned  above. 

For  calculations  in  a,  b,  and  c,  Equations  21  and  22  are  directly  useful,  since  IBM  7094 
programs  have  been  written  for  their  evaluation. 

Additional  thermal  conductivity  measurements  on  tantalum  should  be  made  to  higher 
temperatures  so  that  comparisons  can  be  made  with  the  higher  temperature  literature  data  of 
Figure  12.  Concurrent  electrical  resistivity  measurements  would  also  be  necessary. 

2.  Thermal  Contact  Resistance  Measurement  and  Analysis  Program 

To  allow  convenient  use  of  the  analysis  techniques  developed  in  this  investigation 
two  extremely  time  consuming  steps  should  be  simplified: 

a.  The  analysis  of  the  surface  roughness  and  waviness  data.  And, 

b.  The  analysis  of  the  contact  element  temperature  distribution  data  required  in 
calculating  the  waviness  and  radiation  thermal  resistances. 

Statistical  analyses  of  the  surface  data  as  suggested  by  Henry  (Reference  33)  is  reasonable 
for  random-lay  finishes.  Study  of  practical,  non-Gausian  distributions  should  be  emphasized. 
The  analyses  of  the  computer-generated  finite  element  temperature  fields  could  be  conducted 
with  the  computer  by  incorporating  an  array  of  logic-comparison  operations  in  picking  the 
proper  parameters  for  evaluating  Equations  30  and  31  and  for  determining  A  Tc. 

Thermal  contact  resistance  measurements  to  higher  mean  interface  temperatures  would 
be  extremely  useful  in  studying  in  more  detail  the  radiative  transport  across  different  types 
of  interfaces.  Because  the  contact  conductance  is  not  an  intrinsic  property  at  high  temperatures 
(i.e.  ,  it  is  heat  flux  and  temperature  gradient  dependent),  the  interfacial  flux,  A  Tr  and 
rm  should  be  varied  systematically.  The  equipment  modifications  recommended  above 
would  facilitate  these  measurements. 
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It  is  recommended  that  future  measurements  be  confined  to  surfaces  with  practical 
finishes.  Idealized  contact  elements  such  as  that  depicted  in  Figure  30  have  been 
experimentally  studied  by  many  investigators,  proving  quite  conclusively  that  the  basic  model 
is  physically  sound.  The  major  problem  today  is  developing  usable  ways  of  incorporating 
these  results  in  predictions  for  real  surfaces.  Studies  on  other  nonmetallic  interfaces  are 
necessary  due  to  the  continual  occurrence  of  such  contacts  in  practice  (Reference  1).  Other 
brittle  materials  such  as  aluminum  oxide,  zirconium  oxide,  zirconium  diboride,  hafnium 
carbide,  and  various  silicides  should  be  studied  as  interface  components  along  with  semi- 
brittle  materials  such  as  Union  Carbide  ATJ  graphite  and  the  Poco  Graphite  Corporation 
series  of  materials. 
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APPENDIXES 
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APPENDIX  I 
HEAT  LOSS  ANALYSIS 


The  boundary  conditions  given  in  Figure  8  were  used  to  formulate  the  Laplace  equation 
for  obtaining  the  temperature  distribution  in  the  cylindrical  thermal  insulation  segment  and 
is  expressed  in  the  following  way  for  this  two-dimensional  axisymmetric  case 

V2I  =  +  _£l_  s  o 

d  r 2  r  dr  dz  2 

Implicit  in  the  writing  of  this  equation  are  temperature  independent  properties  for  the  medium 
in  which  the  temperature  distribution  is  sought.  Since  the  differential  equation  is  linear,  the 
influence  of  each  boundary  condition  can  be  determined  separately, that  is,  superposition  can 
be  applied.  The  solution  considering  only  boundary  condition  4  (Figure  8)  is  given  first, 
followed  by  the  solution  considering  boundary  conditions  1  and  2  simultaneously. 

A.  SOLUTION  CONSIDERING  BOUNDARY  CONDITION  4 

The  simplified  boundary  conditions  used  in  this  development  are  as  follows: 


r  -  'o  .  z  .  T  =  0 


r  =  r,  ,  z  ,  T  =  0 


ro  <  t  <  r,  ,  z  =  0  ,  T  :0 

ro  <  r  <  f,  .  z  =  L  .  T  =  f2  (r  ) 


Ts  fz  (r) 
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Assuming  a  solution  of  the  form 

T  (r,  z  )  =  R  <r  )  Z  ( z  ) 
and  substituting  in  the  differential  equation: 


V  2  T 


d*T 


+  -L 

r 


a  T 


a2T 


dr2  r  dr  dz2 

yields  the  following  system  of  ordinary  differential  equations: 


=  0 


__l _  d2  Z(  z  )  jj z 

Z  <  z  >  dz2  =  ^ 


and 


I  d  R(r ) 


R(r)  dr 


I  __ _ I _  dR  ( r  )  w  2 

r  '  R  ( r  >  '  R  (r  ) 


or 


d z  Z  (z  )  ,2 

- -  -  Z  (z  )  K  =  0 

,2  v 

d  z 


and 


d  R  ( r  ) 
dr2 


I  d  R  ( r  )  i/2 

+  /(  R(r)  =  0 


dr 


/  2 


A  question  may  arise  concerning  the  sign  of  \  .  To  satisfy  boundary  condition  4, 

->  necessary  that  the  Fourier  coefficient  integrals  involve  f2  (  r  ).  This  condition  is  met 
when  the  generalized  solution  is  in  the  form  of  the  hyperbolic  function,  zero  order  Bessel 
function  product;  the  sign  of  ^  2  assumed  here  will  lead  to  a  solution  of  this  form.  In  the 
next  portion  of  the  solution  the  opposite  sign  is  required. 

Transforming  the  independent  variables  to  and  i  ,  we  get: 


d  R  (  K  r  ) 


I  dRIX'r  ) 


d(>(r)  d(J(r) 


4-  R  (  f(r  )  =  0 


and 


±J1K£L-  z.*„=o 


98 


J 


AFML-TR-69-1 

The  equation  in  is  in  the  form  of  Bessel's  equation  of  zero  order  and  has  the  general 
solution 

R(4fr)sAJ0(«fr)+BY0(/{r) 

wiiere  J0  and  Y0  are  the  Bessel  functions  of  first  and  second  kind  of  zero  order. 

The  equation  in  has  a  solution  in  terms  of  exponentials  with  real  argument  which, 
in  the  most  general  case,  is  exoressed  as  follows 

Z  ( 1 )  s  c  cosh  +  D  #,nh 

Thus,  a  solution  of  the  original  differential  equation  becomes 

T  (r,z  )  *  [a  J0  (^r  )+BY0  (^r)]  [  Ccosh  (f{i)  +  Q  sinh  (£z)] 

Turning  to  the  boundary  conditions  (Figure  8),  from  B.  C.  3 

- ►  C  =  0 

or 

T  » r,  z  )  :[AJ0  (/(r)  +  BY0  lf(  r  )]  [o  jinh  { z  )  ] 

To  apply  B.  C.  's  1  and  2  the  following  conditions  must  be  met 

[a  J0  (Y'o  >  ♦  BY0  (/(  r„  )  ]  [  0  *.nh  (<z  )  ]  »  f, 

and 

[AJ0  (*r,)  +  BYC  )][  D  s.nh  z  )]  =0 

By  combining  these  equations  and  eliminating  the  constants  A  and  B  the  following  condition 
is  obtained- 

Yo  ‘t'o  1  Jo  {>('<  1  -  Vo  1  ■'o  {K'o  »  =0 

This  condition  is  fulfilled  by  an  infinite  number  of  eigenvalues,  ,  which  are  the 
positive  roots  of  this  equation.  Introducing  cylindrical  function  notation, 

uo(^n-  ro  •  ro>  Yo  <<  r.)  "  V  r.  >  Vfn'o* 

we  can  define  the  eigenvalues 

uo  (^n  .  ro  •  ri  )  =  0  •  n  =  I.  2,3,  . 


l 
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Thus,  there  are  an  infinite  number  of  general  solutions  to  this  boundary  \alue  problem 
of  the  form 

T(r,z  )  =  [AnJ0(^nD  4-  Bn  Y0  r)]  [  Dn  sinh  (^z  )] 

or 

T  lf,z  )  =  [  En  Je  (/fn  r  H-Fn  Y0  r )]  [  s  jnh  (£nz  )] 

Subsequent  Fourier  coefficient  evaluations  are  facilitated  if  the  above  equation  is  cast 

in  terms  of  the  cylindrical  function  (T  ;  r,,r.  ),  i  and  j  =  0  ,1 . 

v  In  i  ) 

Defining  ^  (  Yn  r  )  os 

'  Kn  '  >E  En  Jo  'K„'  I  F„  \ 

then 

<3  1  1=  En  J0  <^n,o  >  +  fo  y0  (*„'„> 

and 

>sEn  J„  <*„'■  1  ■*  Fn  Vo  <*„ 

Combining  these  equations  to  obtain  expressions  for  En  and  Fn  results  in: 

-  f,  >  Yo<*nfo  >  -  [ *n  fo>  Yo  >  _ 

E°  ’  Jo  <Kn',>Yo  <Hnro>-  Jo^n'o  >  V*nf,  > 

fo  >  Yo  (  f'  1 

uo  (*n  '•  ro  •  f.  1 

and 

.  fo  >  ^o  fi>  Jo  <  Yn  re) 

Yo  <Oo>Jo  ,-Yo<\nr.,Jo^nro» 

^  (  'VVi  f  o  )  ^  Vnfi  '  —  ^  ^  n  f  i  ^  Jq  n  ro | 

u0  Kn  >ro  •  r.  > 
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Now, 


)* 


^  K  n  •  ro  »  ri  ^ 


[<?<«'„  'ol-'o  <4fn  !  -G  <*«  I  V^n  '»)]  *„<*„'> 


uo  ^  Kt\  •  ro  »  ri  ) 


_  ^(Vn  r,  >[  Wnfr^o  <  *n  '  >  ~  >*0  (  W  ^0  Ifa  r  )] 
uo  «  ro  *  ri  ) 

f  <g(V„fe>[Jo  <Vn  f.>Yo(*n  r  >-Yo<*r/,>Jo<</  j] 
uo  <  Vn-ro  ’ 


Simplifying,  wo  got 


f>  s 


<^(^n  ri  >[uo(  Vn’ro  ,f  }]+  ^[“o  (^n-r  *r,  >] 


uo  •  ro  >  ri  ^ 


Thus,  a  solution  of  the  differential  equation  can  be  written 

T  (  r,  :  )  =  (%(  r  )  •  smh  ( *  ) 

In  addition,  any  linear  combination  of  solutions  of  the  above  form  will  be  a  solution  of 
the  equation.  That  is, 

T  ( r  ,  z  )  =  Z  (  ^nr>’  ,mh  (^n  z) 

n  - 1 


or 


T  (r,  z)  =  Z 
n  - 1 


>  fo  >]  *  (^>(/lt/nro)[uo(  *_r  >  >f  >]_  . 
uo  <Vn  •  fo  -S’ 


sin h  (/^n  i  ) 


It  now  remains  to  satisfy  B.  C.  4.  By  substitution  in  the  above  equation 


T  (r.L  )  =  f,  (r)=Z 


f,  )  [u0(Vn,r0  ,0]-*-  ^nr0)[u0(Vn  f  J 


n  =  I 


uoKn  .  fo  > 


smh  (^n  L) 
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The  quantities  r, )  and  rQ  )  ,  which  contain  the  constants  En  and  En 

can  be  considered  the  coefficients  of  a  cylindrical  function  Fourier  series  with  a  general 
term  of  the  form  uq  (^  >  r  ;  r,  )  •  smh  (  z )  The  Euler  formulas  will  then  be  of  the  form 


a 


n 


r  f2(  r)  uQ  ( r )  dr 


where  the  norm  N  |  u0  (  r  )  j  is  found  by  applying  Lommel  integral  relations. 


The  final  result  for  the  case  considered  here  is  given  by  Carslaw  and  Jaeger 
(Reference  34)  as: 

T  ( r,  r )  =  ~  £  A.Jii^..fo)uo  ( Yn  ,  ^  r,  Uinh  «nO  j"  f  f?  (f)  Uq(  ^  ,  r ,  r  )  dr 

0:1  [Jo^nro)-Jo<,^nr.,]SinMVnL)  '='o 


where 


and 


(r)  is  on  arbitrary  temperature  distribution  along 
ro  <  r  <  r,  .  L 

are  the  positive  roots  of  uQ  (^n  ,  rQ  ,  r ,  )  =0 


B.  SOLUTION  CONSIDERING  BOUNDARY  CONDITIONS  1  AND  2 

With  the  above  solution  for  boundary  condition  4  now  established,  it  is  necessary  to 
calculate  the  influence  of  boundary  conditions  1  and  2  (Figure  R)  on  the  overall  temperature 
field.  This  is  done  in  a  manner  somewhat  analogus  to  the  above  derivation  except  that  the 
solution  will  be  in  terms  of  modified  Bessel  functions  a^d  trigonometric  functions  to 
accommodate  the  cylindrical  surfaces  boundary  conditions. 
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The  problem  to  be  solved  is  as  follows.  Given,  are  simplified  Boundary  Conditions 
(Specifically,  originally  P.C.'s  1  and  2): 


We  will  assume  a  solution  of  the  form: 

T  (r ,  2  )  -  p  (r  )  ^  iz) 
Substituting  in  the  differential  equation 


(r  ) 


and  rearranging  yield 


dr  2 


+p  ( r  ) 


d y  iz) 
dzZ 


-  0 


_( d2iP(r)  j_  I  dpi r  )  _  _  _j_  dZ  iz  ) 

Pi')  d'Z  +  f  PlO  d  r  '  %iz)  di7 

Since  variables  have  been  separated,  two  equations  result' 


I  d‘P(r) 


Pi  r)  dr 


T 


dr 


'  Pi')  d 
d  2  y  i z  )_  _  _  xz 


i  13  ^  (r).  -  +  x2 


A  =  constont 
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or 


d  V<'l  I 
~~a7r  +  7 


dp(r) 

dr 


-X2  p(t )  -0 


and 


d2  }  CD 
J  * 


(z)  =  0 


\2  .  » 

Again,  a  question  may  arise  concerning  the  sign  of  X  .  If  A"  were  negative,  a 
solution  could  be  obtained  in  terms  of  Jo<r)  and  Y0  l  r)  *  zero  order  Bessel  functions  of  the 
first  and  ?econd  kind,  respectively,  and  in  terms  of  the  hyperbolic  functions  sinh  (  z  )  and 
cosh  (  z  )  ,  as  shown  earlier.  As  indicated  below,  if  X2  is  assumed  positive  the  solution 
is  in  terms  of  IG  ( r )  and  K0(r)  ,  zero  modified  Bessel  functions  of  the  first  and  second 
kind,  respectively,  and  in  terms  of  the  simple  transcendental  functions  sin  (  z)  and  ( i)  cos 
The  latter  alternative  is  used  in  order  that  a  solution  satisfying  boundary  conditions  1  and  2 
can  be  obtained  in  Fourier  series  form. 


The  following  relationships  between  the  alternate  solutions  follows  directly  from 
examination  of  the  differential  equation- 


I0(f  )  =  (.)*  JQ  <  .r ) 


si  n  ( z  )  =  —  s  i  n  h  (  i  z  ) 


-  I 


Transforming  the  independent  variables  to  Xr  and  X  z 

d  2/3  (  X  r  )  l  dp(Xr) 

- r -  -  +  - - -  -  p  (  X  r  )  x  0 

d  (  X  r  )  Xr  d  X  r 

This  is  the  modified  Bessel  equation  of  zero  order  and  has  the  solution: 

P  (X  r  )  =  A'  I0  (  Xr)  4-  b'  K  (  Xr  ) 


and 


d*  y  x 
+  =0 
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This  transcendental  equation  has  a  solution  of  the  form 


|  (\l  )  :  C  coi(Xl)  +  D  *in  (  X  * ) 


Thus,  a  solution  of  the  differential  equation  becomes 


T(f,i):[A'l0(Xr)f  b’ko  (Xr)  J^C'coslXll  +  D^sintXzlj 


From  B.  C.  3- 


C3  *  0 


From  B.  C\  4 


T  (r.z  )  :  [  A*  I0  (Xr)  4-  B*  K0  (Xr  )  J  o'  sm  (Xz) 

T  (r,L)  *  [  A*  X0(Xr)  +  B*  K0  (  Xr  )]  o'  s.n  (  Xl)  -•  0 


This  condition  is  fulfilled  when 


\L  :  rwr,  n  :  o,  I  ■> 


Hence, 


T.,.,|.[E:  i0(-^-)  +  f;k0  (ill-)] 


where  now  the  constants  K1  and  F'  are  subscripted  as  K'  and  F*  ,  since  there  is  a  solution 

n  n 

for  each  value  of  n. 


A  representation  for  Hp  and  Fp  more  convenient  for  Fourier  coefficient  evaluation  can 
be  obtained  bv  defining  the  following  quantities  in  evlindrical  function  form.  In  general,  let 


/»‘Xr>.p(-SfH.  I0  +Fn'K0(^r  ) 


Specifically, 


f  r;  1.(^1^;  K0(^) 
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and 


/  n  rr  r, 

Hi.  v 

(  n  7T  r, 

W  F'  K  / 

f  nv  r,  \ 

V  L 

l  L 

/  +  Fn  Kc  \ 

>  L  / 

Solving  these  equations  for  the  constants  and  results  in 


and 


/  nirri  \  „  /  n  Tr  ro  x 

-  P 

(  n  ^  r0  \ 

*o 

/nTrr,  ^ 

\  l  *  Ko\  L  ' 

'  L  1 

V  L  ) 

l  nwr  ■  \  „  (  n  "’'o  ^ 

T 

1  n  rr  r0  . 

\  K 

l  n7rr.  ^ 

v  l  /  Ko  \  j_  / 

1  0 

1  L 

1  L  > 

n  (-n7rr?')  T  ("*'•  ' 

)  -  o  i 

f  nrrr,  x 

T  / 

n rrr0  \ 

^  V  L  /  0  \  L  - 

/  p  1 

l  L  / 

xo  \ 

.  L  / 

These  equations  can  be  simplified  by  defining  the  function: 

n  v  i.  \  /nr rr.  \  /  n  7r  r„  \  /  n  tt  r, 


I.(^)  K.(^-) 


Thus, 


1 

/nrrr, 

\  K 

fnwr  o\ 

/  n7rr0  \ 

K  ( 

f  n7rri  N 

I 

,  * 
C  -  _ 

V  L 

1  Ko 

\  L  1 

“  P 

\  L  / 

Ko' 

i  L  / 

^  n  - 

Gn<-, 

■  fo  » 

and 

o( 

n7r  ro  \ 

Io( 

n  rr  r,  > 

1  -o 

f  "TT  r, 

)  r  / 

n  tt  r0  ^ 

,  p\ 

r 

L  1 

L  - 

1  P 

l  L 

/  Io\ 

L  1 

* 

Gn 

'  ro 

) 

Finally,  by  writing  the  equation  for  ^  n  ^  f — j  terms  of  these  equations,  the  following 
result  is  obtained- 


/  nrrr  > 

1  P  ( 

,  n  rrr,  ^ 

1  °n  .  r0  >  +  P  1 

f  mrt0  N 

(  L  J 

^  L  / 

P  \  L  1 

G  n  ^  ri  *  r2  ^ 

A  solution  of  the  differential  equation  can  now  be  rewritten  as 

T  ...  (^-) 
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In  addition,  any  linear  combination  of  solutions  of  the  above  form  will  be  a  solution  of  the 
equation.  In  the  most  general  sense, 

T".'|=  I  p(— —)  \~~r) 

n=l 

It  now  remains  to  satisfy  boundary  conditions  1  and  2  along  the  cylindrical  surfaces 
r  =  rQ  and  r  =  r(  .  Appropriate  substitution  in  the  above  solution  yields 

I  p  (^)  „n 
/nfff,  \  /  n  tt  rQ\ 

P  V - \T~)  Gn<V  'o  »  +  P  V— C-9-)  Gn(r. 


=  I 

n  =  I 

Simplifying,  we  get 


G  n  <  r  i  .  ro  ) 


[n-rri  ^ 

'  "n(~T— ) 


n  nr„ 


n  =  l 


(  n  wz 

• 

\  *  y 

/  n7rz  \ 

\  L 

n=l 

an  s,n 

(  L  ) 

r  =  rt: 

(nTri  ) 

■ 

m  y 

£ns,n  ( 

nTT  z  \ 

'  L  / 

*  z. 

n=l 

L  / 

Then,  if  an  and  /3n  are  considered  as  coefficients  of  the  indicated  Fourier  series,  the 
following  expressions  can  be  derived  using  the  standard  procedures  (Reference  2)  in  developing 
Euler  formulas  for  Fourier  coefficients 

«»■«  f  A<.i ..» (-^-) 


and 


a  „  /  nTrr,  \  2  rL  /  n  rr  r  x 

/3„  =  R  (  — - )  *  —  /  I,  ID  sin  (— — )  0, 

Combining  results  obtains 

/  nvz  \ 

_  •  51  n  I - }  .  L 

T  1 1 : 1  ■  r  £ -zXtf  ( G- 1 '  - ^  * 


+  °n  '  1  /  <o  1,1  5,n  ~ )  dl  } 
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This  equation  combined  with  the  result  obtained  separately  for  the  original  boundary  condition  4 
gives  the  following  equation  describing  the  temperature  distribution  in  the  annular  insulation 
sketched  at  the  outset: 


where 


■r ,  ,  ^n^o^nro*“o^n»r'ri^s'n^^nz  C  '  ,  ,  , 

T(r,z)=TZ-— - - - * - J  r  f2(r)u0(Kn,r,rt)dr 

0  =  1  [VKn'o>  -Jo  <Kn',  >Jsi"h  KnL  f° 

,  2  ^  s,n(J!nr)  r  rL  mttz  ^ 

+  Gn  1  •'  >  /*" '»(!)  «i»  ("V  )*“} 

0 

«o  <Kn.S  .  'j  >  z  Jo  <*n',>Yo  <  *n  ' ,  >  “  Jo  <  Kn  'j  >  Yo  <Kn'i  > 

^  *  positive  roots  of  uq(  ^  ;  rQ;  r^  -  0 

V',  .'o'*  K.(-=^) 

fQ(z), fj(z), fgfr)  =  arbitrary  temperature  distributions  subject  only  to  the  requirement 

that  the  above  integrals  exist. 


C.  EVALUATION  OF  FOURIER  COEFFICIENT  INTEGRALS 

An  explicit  equation  for  T(r,  z)  describing  the  temperature  field  in  the  insulation  with 
the  experimentally  observed  boundary  conditions  is  found  through  evaluation  of  the  integrals 
In  the  above  equation,  and  by  making  appropriate  substitutions  for  the  expressions  f0(z), 
fjfz),  and  f2(r). 

For  the  particular  experimental  arrangement  used  here,  f2<r)  is  not  arbitrary.  Actually 
f2<r)  Is  uniquely  defined  by  the  condition  that  V  T=0  along  the  surface,  z  =  L,  where  the 
temperatures  at  positions  r  -  rQ  and  r  -  rj  are  fQ  (L)  and  f i (L),  respectively.  For  the 
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plane  z  *  L,  the  temperature  distribution  defined  by  the  equation  V2  T  =  O  can  be  found 
via  a  one-dimensional  solution 


Rewriting,  we  get 


Separating  variables  and  integrating  produces 


/ 


constant 


or 


dT  .  _£± 

dr  r 


Integrating  again  yields 

/"'  *  - *  T  *  C,  A  r  +  C, 


Evaluation  of  constants  gives 


r  =  r0  ,  T  =  f0  ( L  I 
r  =  r(  ,  T  =  f,  (  L ) 

f0(L)--C,^,  r0+C2  ,  f,  ( L)  =  C,  A  r,  4  C2 

C2  =  f0  (L)  -  C,  JU  rQ  ,  f(  (L)  t  C,  rf  +  fQ  (L)  -  C(  Jn.  rQ 


f,  (L)  -  f0(L  ) 


Ct  =  f0(L)  -Jn,  r0 


t,  CL)-  f0  (L) 

A  ■£■ 
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Thus, 


T  A  r  -M  ,L,-,'(L|-,°IL)  A,, 

Aft)  -M^)  ‘ 

T  =  („  (L)  +  {  (,  (L)  -  f„  (L)}  [  —  f*~  ] 

Consider  now  the  integral  in  the  first  series  of  the  Equation  for  T(r,z),  which  reads  as 


follows: 


rr' 

J  ff,(r)  u0  lKn,r  ,  f,  )df 
r0 

I  ri 

n  (  r,  \  /  {  fo  +  [fi  *  ~  fo  *L  ^  {t^)  }  r  ’  uo*Kn  •  r  •  ri  ^  dr 

AUo)  r° 

r,  rl 

t  [<o(Ll /  r  u 0<Kn.f  .f,  )  dr  4  [f.‘L)-f0(L)]/ru0(Kn;r  ,  r,  (f  ) drl 

A(^)  'o  ro  Vo/  J 

I  r  f  [  *  *  n  ro  )  ~  ^ K  n  r  i  )] 

[fo(L){2 

4-  If,  IL ) -  4  (L)i  I?  ^  f|  ~  ^°l,<n  f|  l..^V  —  -A.fo(g.-,°IK";»|-J°IK^l)}1 

'■  °  n  »«„'  VKn'o>  tk’j0IK„, 0) 
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The  integrals  in  the  second  series  of  the  equation  for  T(r,  z),  for  which  vie  will  assume 
linear  representations  for  both  fQ(z)  and  f^fz),  is  stated  so  that 
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And  substituting  in  the  expression  for  T(r,  z)  we  will  derive  the  following  nomenclature 


m 

7T  r- 

J  Jo(Kn  ro)uo(Kn  -r  .r0)s,nh  (Kn*  *  1 

MrJ  - 

[Jo  "Vo>  -JoIKn  f,l] sin h (K„  L  ) 

[{Tl„  +  -A  Ic.l}  {j0  fK  nr0>  -  J0  «K„  r,  )} 

+  'VloKWo'-^  -J0(Kn',)^  '«}] 

•  I''17’*) 

To  determine  the  total  heat  loss  into  the  insulation  from  the  specimens  it  is  necessary 
to  evaluate  the  radial  temperature  gradient  at  the  surface  r  -  rQ  and  perform  the  integration 

qloss  =27rroki  /  ("a^- )fir  dz 


where 


k|  ■  insulation  thermal  conductivity  (assumed  temperature  and  position  independent  in 
this  development) 


By  writing  T(r,  z)  in  the  following  abbreviated  form, 

T  *r‘  Z  ®  ’  jl  777  [^n  uo  . r  1  ri  )  sm  h  (  K  n  1 )  ] 

'  ro‘ 

m  (nlrZ) 

+  WV  {G"" ' ^  tT'  -Tl, ] +  G",r" " TiF [T»-T'-o  1  "1 } 


where, 


=  Jo  *Kn  ro  J /[{  Jo  (Kn  ro^  ~Jo  (Knri*  }  51  n  h  <  Kn  L)  ] 

77n=T‘-0^(^)('+^ro,-TLI^r  >[Wo>-V*V,>] 

0  ^  ^  n ro ^  ri  "  (  Kn  r,  )  J*,  r0  j 
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Differentiating  this  equation  with  respect  to  r  leads  to  the  following 


'JaT1"^^  i,  £b Vn 


And  evaluating  the  differentials, 
d 


-£-K"V.'  Kn[Yo  «Kn',l  ^  W>- W,1^”  W'] 

Since  ^  is  a  root  of  uq  ( ^  ,  rQ  ‘  r(  1^0 

ro  )  Y0(Knr0)  _ 


^o^nri'  ^  *  n  ri  ^ 


then, 


dr 


Mh.  >]-  4  ^  ^  W  >  -  l  Yo ,v ,] 


Evaluating  at  r  =  r  yields 
o 


d 

dr 


.  =4r  [VKn^>  J>nV-Jo<Vo>r<>,,'0>] 

r-  rn  P 


«n 

r  2 

. _ 2_ 

*n  ri  l  I 

0 

KnroF 

**o  ^  n  ro^ 

For  differentials  of  terms  of  the  second  series  in  T(r,  z), 


_d_ 

dr 


Which  are  defined  as 


then 


o 

M 

M 

n  irr  ' 

IK1 

f  nw'1) 

-ro( 

L  - 

\  L  } 

L 

/  L 

A,  M  L  ) \ 

-  n7r  li 

/  n  ir  r 

k  / 

’  n  x  r0  \ 

+  I0 

/— 

•n  r0 

L  t1' 

\  L 

)  Ko\ 

i  L  / 

V 

L 

)  •  V  L  )/ 

ro  )£  *1  ( 

n  tt  r  i 

)  Ko( 

n  7rr0  v 

+  Ic 

r 

1  K  (  nWf  ) 

L  - 

L  / 

>v 

L 

I  K,\  L  1 

d 

Gn  <' 

•ro  1 

n  tt 

t 

r 

/  _ 

ro’ 

d  r 

L 

Gn 

i f  • 
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Similarly, 


n  it 


Gn  (r.  .  r  »  8  6n  (r  .  f,  > 


dr  " 

where, 

.  /n7rr\  /  n  7r  r.  \  /  n  tt  r  \  /  n  tt  r  \ 

>* x.(— — )  (— r~)  +  «,(— r- ) 

Combining  results  gives 

1  JrI,o  Jl(r-)  n;|  ^  '  •  1  Jo  <Kn'o>  ‘ 

+ -f-  i  { [T  -t  [T  T  (^) 

L  n=|l  Gn(r.  ,ra)  L  '  L!  J  G„(r  r  )  L  °  L0  '  L  / 


Integrating  with  respect  to  Z  to  determine  qjosg  producer 


2 

v  e  _  Jo(^nri! 

f  cos  h  (Kn  L  )  -1  ] 

Mh) 

n*-l  0  "  J0  <Vo> 

l  K,  J 

Alternate  terms  of  the  first  summation  are  zero  due  to  the  factor  [  •  —  (  —  1  )  j  . 
Thus,  only  odd  values  of  n  need  be  considered  in  this  summation. 

In  the  second  summation,  determination  of  tne  eigenvalues,  ^ ^  ,  is  necessary.  These 
values  were  defined  by  the  condition 

^  n  ro  ^  ^  ^  n  ^  —  *^o  ^  ^  n  *  i  )  (  K  n  r© )  :  0 

For  the  experimental  arrangement  used  here, 

-r-~  -  6 
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thus  the  above  equation  can  be  rewritten  in  terms  of  rQ  to  read 


Jo  <Knro)  Yo  16  Kn  ro  >  ”  Jo  <6  Kn  fo  >  Yo  < Kn  ro  >  '  0 


Roots  of  this  equation  for  low  values  of  rj/r0  have  been  tabulated  (Reference  35).  For  the 
case  considered  1  3**e,  the  following  expansion  may  be  used  in  evaluating  the  n^1  root 
(Reference  36) 


'o*n 


2  3 

q-  p  r  -4pq  4-  2  p 


where 


8  = 


n  rr 
~p~^\ 


P 


I 


Bp 


100  (  p*  -  I  )  (32)  (1073)  (/35—  I) 

3(8p)5  ip-  \  )  ’  '  5  (  Bp  )5  {p  -  I  ) 


The  results  are  given  in  Table  I. 


Finally,  the  evaluation  of  qiOS8  requires  the  experimental  temperature  distribution  data. 
This  data  obtained  directly  from  analysis  of  measured  temperature  distributions  is  as  follows: 


TABLE  XI 


BeO  THERMAL  CONDUCTIVITY  GUARDING  DATA 


Data 

Point 

tl 

L1 

T1 

TI 

lio 

T 

o 

Ti 

ki 

2 

221.3 

112.3 

165.  2 

82.5 

166.8 

«■ 

4 

293.0 

131.5 

249.0 

111.0 

212.3 

I3EH 

10 

687.0 

359.5 

748.  0 

165.0 

523.3 

mm  , 

11 

819.3 

421.2 

936.  0 

177.0 

620.3 

ESI  i 

12 

894.2 

458.8 

950.0 

188.0 

676.5 

13 

988.2 

494.8 

972.0 

194.0 

741.6 

0.26 

14 

1060.0 

508.  5 

1215.0 

188.0 

784.3 

0.  29 

16 

1204.2 

533.2 

1537.0 

114.0 

863.2 

0.  ^3 

IS 

1333.8 

541.2 

1910.0 

60.0 

937.  5 

0.  38 
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Physically,  the  surface  r0  <  r  <  r,  ,  i  =0  can  be  assumed  to  be  at  the  temperature 
fQ(0)  =  T  .  Since  in  the  derivation,  it  was  assumed  that  this  surface  was  at  zero  (original 
B.  C.  3),  the  temperature  in  the  above  table  must  be  adjusted  to  this  base. 

The  adjustment  has  no  effect  on  the  calculation  so  long  as  the  val'  ^  of  TTj  are  consistently 
retained. 


TABLE  XII 

ADJUSTED  TEMPERATURE  DISTRIBUTION  DATA 


Data 

Point 

tl 

L1 

T1 

tl 

0 

To 

ki 

2 

138.0 

29.  8 

82.7 

0.0 

0.  08 

4 

172.0 

20.  5 

138.0 

0.  0 

0.  10 

10 

522.0 

194.  5 

583.0 

0.  0 

0.  18 

11 

642.3 

244  .2 

759.0 

0.0 

0.  21 

12 

706.2 

270.  8 

762.0 

0.0 

0.  23 

13 

794.3 

100.8 

778.0 

0.  0 

0.  26 

14 

872.0 

320.  5 

1027.0 

0.0 

0.  29 

16 

1090.2 

419.  2 

1423.0 

0.0 

0.33 

19 

1273.8 

481.  2 

_ 

1850.0 

0.  0 

0.38 
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The  final  form  of  the  equation  for  computer  solution  is  then: 


=  2Tr,0r,[4(T, +TL  )  i  )  £  Gn  1 '°  •  '■ 1 

^  1  nil  nGn  (r,  ,r0)  0  n=  I  n  Gn  ( r,  ,  rQ  ) 


n  =  1,3,5, 


— ~TT  i  T, - - {[TL  (H/..0i-TL  A,o]kIKr,'o'-V<r,', 


*'TL|-TLo'[vv))i',-W,ii-J)vv,'( 


COSh  (K.L)-I 


n  *  I,  2.3. 
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APPENDIX  n 

FINITE  ELEMENT  ANALYSIS 

A.  NODAL  EQUATION  DERIVATIONS 

The  fjllowir^  paragraphs  give  the  derivation  of  equations  for  each  nodal  type  depicted 
in  Figure  32;  the  nomenclature  used  follows  that  given  in  this  figure. 

1.  Internal  Nodes: 

The  internal  nodes  car.  be  represented  as  sketched  at  the  right. 


The  equations  for  the  heat  flux  in  each  direction  are  written 

"If,  om  *  £7  [f  <1  .  J+  II  -  T  (I.J  )]ax 
I,  J  +  I 


q)from  :  ~  [  T  (  I  ,  J  ~  •  1  -  T  (  I  ,  J  )  1  x 

I.J-I  Ay  1  1 


q)from  =  [T<  I-M.Jl-T  (  I,  J  )1  Ay 

X+I.J 


^)from  *  [t  (  I  “I,  J  h-  T  ( I  ,  J  )  1  A  Y 

I- I,  J  a* 
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Summing,  we  get 

[t(  I.  J+l )  +  T  ( I.J-D-2  (I  .J  )]  +  k  -—-[t(1+I.  J  )  +T(I-I.  J)-2  T  (I,  J  )] 


If 


Ax  =  Ay  *  2 


T  ( 1 ,  J  f  I )  +■  T  (I,  J  -  I  )  +  T(  X  +  I  ,  J)  +  T  (l-l)~4[T(ItJ)]  -  O 

2.  Nodes  Along  Lines  of  Temperature  Symmetry 

The  centerline  of  the  conducting  element  (I  *  1,  J)  was  assumed  to  be  a  line  of 
temperature  symmetry  and,  thus,  it  can  be  assumed  that, 

T  ( I +  I , J  )  =  T(l-I.J) 

From  the  result  for  the  previous  case  it  can  be  seen  by  inspection  that 

T  ( 1 ,  J  4-1 )  4-  T  (I,  J-l)  +  T(  X  +  I  ,  J  )  +  T  ( I  -  I  )  -  4  [t  (I,J)]  *0 

or 

T(I .  J )  =  Jj[t<I,J+I)4T(I,J-I  )  -F  2T  C I  ♦  I ,  J  )] 

The  surface  I  ■  9,  J  >  2,  is  an  adiabatic  boundary  and  thus  can  also  be  assumed  to  be  a 
line  of  temperature  symmetry.  In  this  case  the  nodal  equation  will  be  of  the  form 

T(I,J  )=  -^-[T  (I,  J  +  D+T  (I.j-I)  +2T  (I-i.  j)] 

3.  External  Corner  Node 

From  Figure  32  we  see  that  the  external  corner  node  occurs  at  the  position 

I  *  9,  J  =  2.  An  expanded  view  of  this  node  and  definition  of  appropriate  nomenclature  are 

given  below.  The  radiation  heat  transfer  coefficient,  h  is  defined  by 

R 

Qrodiat ion 

hR  At 
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An  explicit  relationship  for  hR  will  be  given  later.  By  writing  the  heat  flux  equations  about 
this  node,  we  obtain 

I.J-I 


<0  =  rb  [T<i.J+l>  - T  (i.j  >]  ax 

1,J+I  AY 


q)  from 

r  +  i.  J 

k 

Ax 

L  T<I  +  I.J  ) 

-  T  (I  . J) 

AY 

2 

^  )  f  rom 
I-I.J 

k 

T(1  -1,  J  ) 

—  T  (  I ,  J  ) 

AY 

Ax 

2 

Since  the  surface  I  s  3,  J,  is  insulated, 

T  (  I  —  I,  J  )  =  T  (I  +-  I  ,  J  ) 

Summing  the  heat  fluxes  yields 

hR  AX  [m,  J-l  )  -T  <1,J)]  +  -*  *-  [  T(I.J+  I)  -  T  ( I.  J  )] 

+  2  ( ~Ex  '  <  i -i.J)  -t(  r.  j )]  =o 
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Simplifying,  we  get 


(  T  ( I  -I .  J)AY2+T(I,  J-H)AXZ  +  4s-  Ax2Ay  T(I,J-I)  > 

t(i,j  )  =  \ - 7—r-. - f 

1  r  -  Ax 2  Ay  1  ' 


[ax?  ■♦•Ay5 


Since 


Ax  =  Ay  =  £  then 


T(I  -I.  Jlf  T  (  I  .J-M  )  + 


hR  £ 


T  (I, J-l ) 


T  (  I  ,J  )  -- 


4.  Internal  Corner  Node 

Referring  again  to  Figure  32  we  see  that  an  internal  corner  node  occurs  at  I  ■  4, 
J  -  2.  Definition  of  terms  for  this  node  are  below. 


The  heat  flux  equations  about  this  node  are  written  as  follows 

»  Y  [  T  (r  ~  1.  J  1  ~  T  ( I  ,  J )]  AY 
t-I,J  a  J 
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")f,om  =  3^- [  t*i  +  ■.  j>  t  ii.ji]  ~ 


I  +  i,  J 


«)i,om  =  -Sy  [tu.J+h-  t  ii.ji]  ax 


T.  J  +  l 


q )  ( 10m 

I.J-I 


hR  [  T(  I.  J  -  I  )  -  T  (  L,  J  )]  Ap- 
-h  —  [t  (  I.J-I  )  -T  (  I.J)  ] 

Av  L  1  7 


Summing  produces 


[r|I-l,J|-TII,Jl]*||j  [t(IH,JI-T(I,J)]+5|^-[t(I.J«)-T(I.J|] 


h0  AX 


[ki.  j-i)-t(i,j)]  +  -  0 


Simplifying  gives 


r  i?r  i  ,  hRAx  AY 

,  [?T(I-l,J)+T(I+l,J)jAY  ♦[2T(r,J+l)  +  T(I,J-l)]AX  ♦ - ---T(I.J-I), 

3^Ax2  +  Ay2  J  ♦  f  Ax2  Ay  j 


And,  since 


Ax  =  Ay  -  £  j  then 


T  (I  ,J  )  r 


r  hp  x 

2[T(I-|,J)f  T(I,  J+I1+  T(I  +  I,J)  +  T(I  ,J-I)+  — - -  (I.J-D 

h  o 


5.  Radiation  Nodes 


Radiation  nodes  occur  at  positions,  I  *  5-8,  J  *  2,  and  are  characterized  by 
radiation  heat  transfer  only  in  the  v-direction  from  the  surface  J  =  2  to  the  constant 
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temperature  surface  J  -  1.  A  sketch  of  the  heat  transfer  paths  is  given  below. 


'T  (I.  Jfl) 


T  CX-I.J  » 


Tll-fl.  J)  'Z 


T(I.J) 


T(L.J-I) 


The  flux  equations  are: 


1 )  from 
I.  J4I 

*  )  from 

I.J-I 


I-I.J 

^Jfrom 

141, J 


=  X7  [t(i.J4  I )  —  T  C  J,  J  ]  AX 

=  hR  [t  { r .  J  -I  )  —  T  (  I , J  )]ax 
k  r  i  Ay 


Summing,  we  get 


W'1- 


J4l)-T(r,J)]4^^-[T(r-l,J)— T(I4I,J  )-2T(  I,  J)] 
4  hR  A  X  [t(  I  ,  J  -I)  -  T  (  I ,  J  )  ]  =0 


T  ( I .  J  )  = 


r  1  2  .  2hR  Ax  AY 

[T(I-l,J)4T(l4l,J)jAY  4  2TU.J4I  )AX  4 - - -  T(I,  J  —  I) 

[2  (  AX  2  4  AY2  )  4  AY  ] 
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Since 


AX  =  AY  =  J. 


T  (I,J  )  = 


2hR  1 

T(I’l,J)fT  (I  +  I,J)  +  2T{I,J+I)  f  ^ -  T  ( I ,  J  —  I ) 

__ 


6.  Heat  Transfer  Coefficient  Formulation 


From  the  definition  of  hR  as  qR  /AT  ,  it  follows  that, 


(  T.  -T„  ) 


where 


=  geometrical  view  factor 
=  emittance  view  factor 
*  temperatures  of  radiation  surfaces, 


Setting  T,  :  T0  +  AT  we  then  have 


[T,4  ~To  ]=  [(To  +An  -To  ] 


=  4tq3  At  4  6T02  at  +4Tb  At3  +-At4 


We  can  write 


A  T  •  (  T,  -  T0  ) 


and  simplify  the  above  to 


[V  -t;j 


T.5  +  To  T,  To  J‘  +  T,’ 


Thus, 


+T!2  ro  +T,  To  +Tol 


z  .  3 


'R  '  AT 


Table  VI  gives  a  summary  of  the  necessary  nodal  equations  derived  above. 


124 


1 


AFML-TR-69-1 

B.  NUMERICAL  EVALUATION 

The  table  below  gives  a  "first-guess"  temperature  distribution  based  on  the  flux  pattern 
depicted  in  Figure  31,  assuming  that  the  surfaces  I,  J  =  1,  and  I,  J  =  8,  are  fixed 
temperature  planes  at  500  and  1500*  C,  respectively. 


TABLE  XIII 

"FIRST  GUESS"  TEMPERATURE  DISTRIBUTION 


Node 

8 

1500 

1500 

1500 

1500 

1500 

1500 

1500 

1500 

1500 

7 

1400 

1400 

1400 

1400 

1400 

MOO 

1400 

1400 

1400 

6 

1300 

1300 

1300 

1300 

1300 

1300 

1300 

1300 

1300 

5 

1200 

1200 

1200 

1200 

1200 

1200 

1200 

1200 

1200 

4 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

3 

880 

890 

900 

920 

950 

980 

1010 

1040 

1070 

2 

700 

710 

720 

740 

790 

840 

900 

960 

1020 

1 

500 

500 

500 

500 

500 

500 

500 

500 

500 

J/I 

1 

2 

3 

4 

5 

6 

7 

8 

9Node 

Assigning  a  quantitative  fceometrical  size  to  the  nodal  grid  (i.e. ,  defining  the  magnitude 
of  J!  in  the  case  where  A  X  :  AY)  is  dependent  on  the  magnitude  of  the  heat  flux  expected 
for  the  assumed  overall  temperature  drsp: 

AT  =  1500  -  500  =  1000°  F. 

Our  experimental  data  has  shown  that  a  heat  flux  of  q/A  =  30,000  BTU/hr-ft^  is  a  typical 
value.  If  we  assume  a  mean  conductivity  value  of  400  BTU-in/hr-ft?  -*  F,  which  was 
observed  experimentally  for  BeO  at  1000°  F,  the  following  equation  results, 


AY  =  AX  -1  £ 


X  A  T  _  (400)  ( 1000  ) 
(q/A)  30,000 


i  S2  0  in 
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f 

Assuming  a  total  hemispherical  emittance  of  t  *  0.  3  for  BeO  will  define  the  following 
emittance  shape  factor: 


y  z  - ! -  =  0  18 

'  [4-'] 

For  this  calculation  it  will  be  assumed  that  -1.0.  Since  T0  -  50(f  F  -  960“  R  and 

Tj  is  a  variable,  the  final  equation  for  hpp  becomes 

hRs  H  T(  )  =  {3  33  XIO'10  T(5  4-  3.20  X  10*7  T2  +  3  07  XI0'4  T(  +0  29a} 

Inserting  numerical  values  for  and  T(I=1,8,J=1)  =  50(7  F  into  the  appropriate 

nodal  equations  leads  to  the  set  of  equations  displayed  in  Table  XIV  for  computer  evaluation 


TABLE  XIV 

NODAL  EQUATIONS  FOR  COMPUTER  EVALUATION 


Type 

Nodal  Equation 

Internal 

Same  as  Table  VI 

Line  of  Symmetry 

Same  as  Table  VI 

Internal  Corner 

External  Corner 

T(I-I,J)  +  T  {  L,  J  +  1 )  +  2  5hR 

T  (  T  .  J  )  =  - — 

2.0  +0.005hR 

Radiation 

..  T(I-I,  J)+  T(I4l,J)+2T(I,.i  +  l)+5.0hR 

T  (  I  ,  J  )  - 

4.0  +0  01  h r 

Heat  Transfer  Coefficient 

h(I,J)=333X|0''°T5f3.20X|0'7TZ  +  3.07X  io"4T  -(-  0294 
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APPENDIX  III 

PHYSICAL  PROPERTY  DATA 

A.  THERMAL  CONDUCTIVITY  OF  THE  ALUMINA  BUBBLE  INSULATION 

To  determine  the  heat  exchange  between  the  specimen  column  and  the  alumina  bubble 
insulation  using  Equation  22,  it  was  necessary  to  have  values  of  the  insulation  conductivity. 

This  conductivity  was  not  measured  in  the  investigation;  thus,  it  was  necessary  to  estimate 
values  from  data  in  the  literature. 

Literature  data  on  alumina  powders  are  reported  for  a  range  of  densities  and  particle 
sizes.  Extrapolation  of  the  data  to  other  ranges  of  density  or  particle  size  is  not  straight¬ 
forward  because  of  the  high  radiative  contribution  to  the  total  thermal  conductivity.  The  data 
of  Laubitz  (Reference  37)  was  found  to  be  particularly  useful  since  the  solid  conduction  and 
radiation  contributionc  to  the  effective  conductivity  were  presented  separately  as  a  function  of 
temperature.  Further,  the  data  were  found  to  be  consistent  with  other  alumina  powder  data 
in  the  same  density  and  particle  size  range  (References  38  and  39). 

The  results  of  Laubitz,  measured  in  vacuum,  on  a  powder  with  density  of  1.77  gms/cm3 
and  particle  sizes  in  the  range  300-420  is  given  as  curve  (F)  in  Figure  41.  This  particle 
size  range  encompasses  the  average  particle  diameter  of  the  alumina  bubbles  used  here 
d  =  396  fj.  .  The  alumina  bubbles  have  a  lower  density,  however,  the  average  value  being 

o 

1.01  gms/cm  .  The  bubbles  are  hollow  but  sii.^c  the  radiative  contribution  is  essentially 
size  dependent  (See  Equation  2  of  Reference  37),  it  was  assumed  that  this  contribution  was 
the  same  as  that  inferred  by  Laubitz  for  material  having  the  same  particle  size  range  as  the 
alumina  bubbles  (curve  (z) ,  Figure  41).  The  solid  conduction  contribution  for  the  alumina 
bubbles  was  calculated  by  reducing  that  measured  by  Laubitz  in  the  ratio  1.01/1.77.  The 
estimated  effective  conductivity  for  the  alumina  bubbles,  found  by  adding  the  conductive  and 
radiative  contributions,  is  given  as  curve  (3)  in  Figure  41.  That  this  curve  is  of  the  right 
order  of  mangitude  is  demonstrated  by  the  dotted  curve  in  the  figure.  This  curve  gives  vacuum 
environment  conductivity  results  by  Wechsler  and  Glaser  (Reference  40)  for  an  alumina 
bubble/graphite  fiber  composite. 

Data  points  taken  from  curve  (3)  are  tabulated  in  Table  XI,  Appendix  I,  where  they 
were  used  in  calculating  the  heat  exchange  between  the  specimen  column  and  the  insulation. 
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Figure  41.  Thermal  Conductivity  of  Aluminum  Oxide  Powder  Insulations 
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B.  BERYLLIUM  OXIDE  EMITTANCE  DATA 

Total  hemispherical  emittance  data  is  required  in  the  evaluation  of  the  emissivity  shape 
factor,  Equation  33.  To  provide  this  data  spectral  normal  emittance  measurements  were  run 
on  beryllium  oxide  disks  using  the  equipment  and  methods  described  in  Reference  30.  A  ratio 
recording  Perkin-Elmer  model  13  U  spectrophotometer  was  used  in  double  beam  mode.  The 
basic  measurement  consisted  of  a  comparison  betwe  n  the  emissive  power  of  the  BeO  specimen 
surface  and  the  emissive  power  of  a  blackbody  reference  cavity  at  the  same  temperature. 
Systems  calibrations  were  performed  with  National  Bureau  of  Standards  reference  standards. 

A  sodium  chloride  prism  was  used  to  obtain  spectral  normal  emittance  data  over  the  range 
from  1.8  to  14  microns.  The  results  at  a  specimen  temperature  of  980°  F  are  shown  in 
Figure  42  for  two  BeO  specimens  of  different  surface  roughness  and  density.  The  specimen 
of  higher  surface  roughness  has  higher  spectral  emittance,  particularly  in  the  region  where 
the  absorption  coefficient  is  high  (4  to  9/i),  and  where  emission  is  predominately  a  surface 
phenomenon.  This  behavior  is  consistent  with  the  expected  phenomenon  in  the  geometrical 
optics  regime:  increased  emittance  with  increased  surface  roughness.  Also,  the  data  is 
typical  of  a  dielectric  material  in  the  infrared  and  compares  favorably  with  other  published 
data  on  the  spectral  emittance  of  beryllium  oxide  (Reference  41). 

The  spectral  data  was  converted  to  total  normal  emittance  form  by  use  of  the  following 
equation, 


* 


TN 


Ebb  (X)  dX 


(X)  dX 


^2 

?  *XN  EBB  ^ 

_Ai _ 

X2 

£  ebb  ^ 


Numerical  integration  was  performed  using  the  spectral  data  in  Figure  42,  =  MX), 

and  blackboly  emissive  power  curves  calculated  from  the  universal  Planckian  distribution 
given  in  Ref<  re  nee  32.  To  obtain  an  value  at  a  given  temperature,  the  blackbodv 

distribution  at  that  temperature  is  used  in  the  above  equation.  Since  the  curves 

were  obtained  at  a  constant  temperature  of  980*  F  it  is  implicitly  assumed  that  the  spectral 
emittance  of  BeO  is  not  strongly  temperature  dependent.  Experimental  evidence  available 
(References  41  and  42)  indicates  that  the  spectral  emittance  is  nearly  temperature  independent 
and  thus  the  evaluation  procedure  is  correct.  The  final  results  of  the  numerical  integrations 
are  presented  in  Figure  43  as  versus  T.  Literature  data  on  similar  BeO  materials 
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Figure  42.  Spectral  Normal  Emittance  of  Beryllium  Oxide  in  the  Infrared 
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is  also  presented  in  this  figure.  As  the  temperature  increases,  the  total  normal  emittance 
would  be  expected  to  decrease,  because  the  blackbody  curve  is  shifting  toward  shorter 
wavelengths  as  indicated  in  Figure  42,  weighing  the  spectral  data  in  the  1  to  4  micron  range 
more  heavily. 

A  small  correction  is  required  to  convert  total  normal  to  total  hemispherical  emittance 
for  this  dielectric  material,  which  has  a  zero  absorption  index.  The  correction  curve  given 
by  Dunkle  (Reference  43)  was  used  to  produce  final  curves.  The  equation  of  the  mean 

curve  used  in  the  computer  evaluation  of  Equation  33  was  as  follows: 

200*  F  <  T  <  700*  F 

«TH  -  0.650  -  (T-200)  (0.000023) 

70 0*  F  <  T  <  2000*  F 

<TH  -  0.639  -  (T-700)  (0.0000129) 

The  characteristics  of  the  BeO  specimens  on  which  total  normal  emittance  literature 
data  is  reported  in  Figure  43  are  summarized  in  Table  XV.  There  is  no  clear  rationale  for 
the  wide  disparity  in  results,  although,  for  dielectric  materials,  both  surface  and  volume 
emission  are  important  over  fairly  narrow  wavelength  intervals.  This  implies  that  the  total 
emittance  will  be  a  strong  function  of  specimen  surface  finish,  chemical  purity,  and  density. 
The  influence  of  density  is  particularly  important  for  near-theoreMcal  density  materials, 
since  the  porosity  (and,  in  turn,  the  effective  optical  scattering  cross  section)  changes 
rapidly  in  this  range.  From  the  characterization  information  available  in  the  literature, 
it  is  clear  that  the  rsults  presented  in  Figure  43  are  on  BeO  materials  of  widely  differing 
character.  Incompleteness  of  the  supporting  data  prevents  quantitative  rationalization  of 
these  differences. 
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Figure  43.  Total  Normal  Emittance  of  Beryllium  Oxide  as  a  Function  of  Temperature 
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APPENDIX  IV 

TABULATED  EXPERIMENTAL  DATA 


A.  BERYLLIUM  OXIDE  THERMAL  CONDUCTIVITY  DATA 

The  thermal  conductivity  data  given  below  was  obtained  bv  comparing  the  temperature 
gradients  in  the  Armco  iron  meter  bars  with  the  temperature  gradient  in  the  beryllium  oxide 
specimen  (Figure  9).  Since  the  conductivity  of  the  Armco  iron  is  known  accurately,  the 
following  equation  can  be  used  to  calculate  the  BeO  conductivity.  The  kupp(lr  and  k]ower 
values  for  the  BeO  were  determined  using  the  temperature  gradient  in  the  upper  and  lower 
Armco  Iron  specimen,  respectively. 

kb,o  *  KAr  [(AT /Ax>a  j  /(AT/Ax:0eO] 

All  conductivity  values  have  the  units  BTU-in. /hr-ft*--"  F.  The  k  column  is  the  arithmetic 
average  of  the  previous  columns. 


TABLE  XVI 


BeO  THERMAL  CONDUCTIVITY  RESULTS 


Data  Point 

Mean  Temperature,*  F 

k 

upi*-r 

**  lower 

O 

1 

97.3 

1584. 0 

1602.0 

1593.0 

2 

129.  1 

1522. 0 

1444.0 

1483. 0 

3 

156.9 

1478. 0 

1397. 0 

1437.5 

4 

181.3 

1353. 0 

1322.0 

1337. 5 

5 

235.2 

1206.0 

1168.0 

1187.0 

6 

296.  2 

1 167. 0 

1106.0 

1136.5 

7 

364.9 

1032.0 

973.  0 

1002. 5 

8 

403.7 

979.0 

930.  0 

954.  5 

9 

501.  1 

818.5 

775.  0 

796.8 

10 

602.  5 

716.5 

673.  5 

695.  0 

11 

618.3 

679.  0 

642.  5 

660.  9 

12 

634.2 

652.  5 

630.0 

64 1 . 3 

13 

758.9 

559.  0 

535.5 

547  .3 

14 

860.9 

512.  5 

471.5 

492.0 

15 

933.  1 

431. 0 

434.0 

432.  5 

16 

979.4 

400.  0 

411.5 

405.8 

17 

1032. 9 

370.  0 

389.0 

379.  5 

18 

1 1 15.  y 

340.8 

360.  2 

350.  5 

IS 

718.7 

520.0 

566.  5 

543  .  3 

20 

562.2 

656  .0 

697.5 

676.8 

21 

565.4 

673.  5 

701.0 

687.3 

22 

248  .  5 

1088. 0 

1098.0 

1093.0 
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B.  TANTALUM  THERMAL  CONDUCTIVITY  DATA 

The  tantalum  data  given  in  Table  XVII  was  obtained  in  the  same  manner  as  that  described 
for  the  beryllium  oxide. 


TABLE  XVII 

TANTALUM  THERMAL  CONDUCTIVITY  RESULTS 


Data  Point 

Mean  Temperature,*  F 

k 

upper 

^lower 

*Ta 

1 

125.6 

371.4 

412.2 

391.8 

2 

203.  2 

386.  5 

424.5 

405.5 

3 

206.  0 

37  1 .  5 

426.3 

398.9 

4 

300.  5 

384  .  0 

421.7 

402.9 

5 

386.6 

391.  2 

415.7 

403.  5 

6 

523.  2 

391.  5 

407.4 

399.  5 

7 

689.4 

385.  5 

392.9 

389.2 

8 

837.  1 

395.  0 

397.3 

396.  2 

9 

1065.4 

382.  1 

396.8 

389.5 

C.  THERMAL  CONTACT  RESISTANCE  DATA 

Tables  XVIII  through  XXVII  give  the  tabular  experimental  data  obtained  on  this  program 
during  the  contact  resistance  measurements.  The  second  column  of  each  table  gives  the 
interfacial  temperature  drop.  The  next  two  columns  give  the  heat  flux  in  the  specimens  above 
and  below  the  interface,  respectively,  while  the  following  two  columns  give  the  corresponding 
thermal  contact  conductance  values.  The  next  column  gives  the  mean  interface  temperature. 
The  last  column  is  the  average  thermal  contact  conductance  calculated  as  the  arithmetic  mean 
of  the  conductance  values  in  the  previous  columns. 
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THERMAL  CONTACT  CONDUCTANCE  RESULTS  -  40  fJ.  in./40^  in 
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THERMAL  CONTACT  CONDUCTANCE  RESULTS  -40  /u  in.  /40 
INTERFACE-  368.6  psi  COMPRESSIVE  LOAD 
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